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Modern  airport  surveillance  radars  (ASR),  situated  on  or  near  most  major  air  terminals,  feature 
coherent  pulse-Doppler  processing,  a  vertical-fan  beam  and  rapid  asimuthal  antenna  scanning  for 
detection  and  tracking  of  aircraft.  Those  radars  might  serve  an  additional  useful  role  by  making 
radial  wind  measurements  in  the  Immediate  vicinity  of  an  airport  so  as  to  provide  data  on 
thunderstorm  outflow  winds.  This  report  presents  a  preliminary  analysis  of  the  capabilities  and 
limitations  of  ASRs  in  measuring  outflow  winds,  Principal  results  are;  (1)  radar  sensitivity  is 
adequate  to  measure  winds  associated  with  weakly  reflecting  (5-20  dBZ)  thunderstorm  outflows  at 
ranges  less  than  20  km  provided  that  appropriate  operating  parameters  are  chosen;  (2)  overhanging 
precipitation,  often  moving  at  a  markedly  different  radial  velocity  than  the  outflow,  will  he  a 
significant  source  of  interference  owing  to  the  vertical-fun  antenna  pattern.  If  radar  reflectivity  is 
approximately  constant  with  altitude,  this  interference  will  limit  the  maximum  range  for  reliable 
outflow  velocity  measurements  to  about  20  km  for  an  outflow  that  extends  1000  m  above  the  surface 
and  to  7  km  for  an  outflow  thut  extends  only  300  in  above  the  surface;  (3)  At  two  example  major  ulr 
terminals  (Memphis  International  and  Denver  Stapleton)  ground  clutter  suppression  of 
approximately  AO  dB,  combined  with  the  use  of  inter-clutter  visibility  techniques,  would  result  In 
an  adequate  slgnal-to-interference  ratio  for  thunderstorm  outflow  velocity  measurement  over  the 
significant  upprouch/depurture  corridors.  This  result  applies  when  the  radar  reflectivity  factor  in  the 
outflow  Is  20  dBZ  or  greater  und  the  associated  winds  extend  at  least  300  nt  above  the  surface. 
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ABSTRACT 


I 


Modern  airport  surveillance  radars  (ASR),  situated  on  or  near  most  major 
air  terminals,  feature  coherent  pulse-Dopplcr  processing,  a  vertical-fan  beam  and 
rapid  ailmuthal  antenna  scanning  for  detection  and  tracking  of  aircraft.  These 
radars  might  serve  an  additional  role  by  making  radial  wind  measurements  in  the 
immediate  vicinity  of  an  airport  so  as  to  provide  data  on  thunderstorm  outflow 
winds  for  use  with  a  low  level  wind  shear  alert  system  (LLWAS)  and/or  a  termi¬ 
nal  Doppler  weather  radar  (TDWR).  The  feasibility  of  radial  wind  measurements 
with  ASRs  is  being  investigated  through  numerical  simulation  and  a  held  meas¬ 
urement  program. 

This  report  presents  a  preliminary  analysis  of  the  capabilities  and  limitations 
of  ASRs  In  measuring  outflow  winds.  Principal  rosults  arc; 

(1)  Radar  sensitivity  is  adequate  to  measure  winds  associated  with  weakly 
reflecting  (5-20  dBZ)  thunderstorm  outflows  at  ranges  less  than  20  km  pro¬ 
vided  that  appropriate  operating  parameters  are  chosen; 

(2)  Overhanging  precipitation,  often  moving  at  a  markedly  different  radial  velo¬ 
city  than  the  outflow,  will  be  a  significant  sourco  of  interference  owing  to  the 
vortical-fan  antenna  pattern.  If  radar  reflectivity  is  approximately  constant 
with  altitude,  this  Interference  will  limit  the  maximum  range  for  outflow 
velocity  measurements  to  about  20  km  for  an  outflow  that  extends  1000  m 
above  the  surface.  For  an  outflow  that  extends  only  300  m  above  the  sur¬ 
face,  reliable  velocity  measurements  will  be  possible  to  about  7  km  in  range. 
Thunderstoi  n  outflows  that  extend  well  above  1000  in,  such  as  gust  fronts, 
could  be  reliably  measured  at  considerably  greater  ranges; 

(3)  At  two  example  major  air  terminals  (Memphis  International  and  Denver  Sta¬ 
pleton)  ground  clutter  suppression  of  approximately  40  dl3,  combined  with 
the  use  of  inter-clutter  visibility  techniques,  would  result  in  an  adequate 
signal-t.o-interfercncc  ratio  for  microburst  or  gust  front  velocity  measure¬ 
ments  over  the  significant  approach /departure  corridors.  This  result  applies 
when  the  radar  reflectivity  factor  in  the  outflow  is  20  dI3Z  or  greater  and  the 
associated  winds  extend  at  least  300  m  above  the  surface. 
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A  Preliminary  Assessment  of  Thunderstorm  Outflow  Wind 
Measurement  with  Airport  Surveillance  Radars 


I.  INTRODUCTION 

The  Federal  Aviation  Agency  is  procuring  over  130  new  airport  surveillance  radars  (ASR-9), 
Their  primary  mission  is  to  detect  aircraft  while  rejecting  ground  and  precipitation  clutter.  A 
separate  weather  channel  will  provide  displays  of  precipitation  reflectivity  for  use  in  air  traffic  con¬ 
trol.  In  view  of  their  high  data  update  rate  (4.8  seconds/scan),  coherent  pulse-Doppler  processing 
and  on  airport  location,  it  is  appropriate  to  consider  whether  airport  surveillance  radars  might  also 
detect  thunderstorm-generated  low  altitude  wind  shear  (LAWS)  such  as  mlcrobursta  or  gust  fronts. 
This  would  require; 

(1)  estimation  of  the  radial  wind  velocity  associated  with  the  outflows; 

(2)  an  adaptation  of  the  LAWS  detection  algorithms  that  are  being  developed  for  use  with 
pencil-beam  Doppler  weather  radars  (1,2,3). 

Figure  1-1  sketches  a  vertical  cross-section  of  a  microburst.  A  microburst  occurs  when  an 
intense,  small  scale  downdraft  encounters  the  earth's  surface,  producing  a  brief  (<  10  minute)  out¬ 
burst  of  highly  divergent  horisontal  winds  (4).  Aircraft  taking  off  or  landing  through  a  micro- 
burst  experience  increased  lift  as  they  enter  the  outflow,  followed  by  a  reduction  In  lift  in  the 
downdraft  and  tailwind  sections  of  the  microburst.  Owing  to  the  small  spatial  extent  of  micro¬ 
bursts  (  <  4  km)  this  sequence  is  encountered  in  a  period  of  1  minute  or  less.  For  aircraft  on 
approach  to  landing,  an  attempt  to  keep  the  plane  on  its  glide  slope  by  reducing  power  in  the 
headwind  section  of  the  microburst  may  compound  subsequent  loss  of  lift  in  the  downdraft  and 
tailwind  sections.  A  r  mber  of  commercial  carrier  accidents  in  recent  years  have  resulted  from 
microburst  wind  shear  (4,5,8). 

Gust  fronts,  as  sketched  in  Figure  1-2,  are  generated  by  thunderstorm  outflows  that  persist  over 
a  longer  period  of  time;  thus  the  leading  edge  of  the  outflow  may  be  distant  from  the  generating 
storm.  Squall  lines  such  as  encountered  in  the  midwestern  portions  of  the  United  States  produce 
gust  fronts  that  propagate  tens  of  kilometers  away  from  the  generating  precipitation.  The  wind 
shear  encountered  by  an  aircraft  penetrating  a  gust  front  is  generally  considered  less  hasardous 
than  that  associated  with  a  microburst  since  the  change  is  toward-  greater  lift.  However,  tur¬ 
bulence  behind  a  gust  front  can  be  hasardous  and  the  change  of  wind  direction  behind  the  front  is 
of  concern  for  runway  usage,  Gust  fronts  are  longer  lasting  phenomena  (15  minutes  or  more)  and 
may  result  in  significant  disruption  of  normal  airport  operation.  Tracking  and  prediction  of  gust 
front  arrivals  at  an  airport  based  on  low  altitude  radial  velocity  data  appears  feasible  (1)  and 
would  yield  Important  cost  benefits  for  airport  operation  (8). 

The  limited  number  of  op  nationally  oriented  LAWS  detection  programs  to  date  indicate  that 
the  regions  of  greatest  practical  concern  for  LAWS  detection  are  one  nmi  wide  corridors  centered 
on  the  runway  centerlines  and  extending  out  about  3  nmi  beyond  the  runway  ends  (8).  These  are 
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the  region*  in  which  aircraft  on  approach  or  takeoff  arc  less  than  1000  feet  above  ground  level, 
LAWS  in  the  area  surrounding  these  corridors  Is  of  concern  because  the  area  of  wind  shear  may 
grow  or  move  into  the  critical  region.  Figure  1-3  sketches  the  key  areas  for  microburst  detection  at 
a  hypothetical  airport.  Detection  of  gust  fronts  at  distances  up  to  25  mni  would  be  useful  for  pro¬ 
viding  predictions  of  significant  wind  shift  at  an  airport. 

The  primary  wind-shear  sensor  under  development  is  a  dedicated,  terminal  Doppler  weather 
radar  (TDWR)  that  will  use  its  pencil  beam  in  a  combination  of  PPI  and  RHI  scans  to  provide 
volumetric  measurements  of  storm  reflectivity  and  velocity  parameters.  The  TDWR  will  provide 
the  capability  for  nowcasting  of  ful’y  developed  LAV/S  .vents  and  ior  identification  of  potential 
microburst  precursors  aloft  as  well  as  detection  and  tracking  of  gust  fronts,  areas  of  heavy  precipi¬ 
tation,  turbulence  and  tornadoes  (9). 


It  has  been  suggested  that  airport  surveillance  radars  could  serve  as  useful  adjunct  LAWS  sen¬ 
sors  for  airports  that  will  receive  a  TDWR  and  as  stsnd-aiono  systems  at  those  terminals  that  are 
not  slated  for  a  TDWR.  The  siting  of  an  ASR,  often  at  runway  Intersections,  or  at  one  end  of  a 
major  runway,  might  allow  for  an  improved  measurement  of  headwind-tailwind  shear  when  a 
TDWR  Is  sited  off-airport.  It  is  unlikely  that  TDWRs  will  be  deployed  extensively  at  secondary 
terminals  or  in  locales  where  LAWS  is  infrequent.  At  these  terminals,  an  ASR  might  provide  at 
least  some  capability  for  wind-shear  detection,  particularly  in  the  critical  region  within  5  ■■  10  km 
of  an  airport.  It  might  also  be  possible  to  improve  the  performance  of  the  low-level  wind  shear 
alert  system  (LLWAS)  by  using  the  data  from  an  airport  surveillance  radar  to  reduce  the  LLWAS 
false-alarm  rate  and  to  distinguish  between  various  forms  of  wind  shear. 

Table  1-1  summarises  principal  parameters  of  the  ASR-9.  The  dual-beamed  ASR-8  currently  in 
operation  at  many  airports  has  very  similar  radar  system  parameters.  If  airport  surveillance 
radars  are  to  be  used  for  LAWS  detection,  the  most  likely  method  would  be  by  way  of  a  processing 
retrofit  that  would  not  affect  their  primary  aircraft  detection/tracking  mission.  Thus  wind-shear 
detection  with  an  ASR  will  be  limited  by  such  factors  as  the  fan-shaped  elevation  antenna  pattern, 
the  short  coherent  processing  intervals  (CPI)  and  the  concomitant  difficulties  In  achieving  adequate 
signal-to-interference  ratios  for  estimation  of  radial  velocity  in  shallow,  possibly  low-reflectivity 
microbursts  or  gust  fronts. 


Table  1-1 

ASR-9  Parameters 

Frequency 

Trantmitter 

2. 7-2. 9  GHs 

Polarisation 

Linear  or  Circular 

Peak  Power 

1.1  MW 

Pulse  Width 

1.0  nn 

Block-Staggered  CPI  lengths 

8  pulses/ 10  pulses 

PRFs  (Example) 

972  * -1  / 1250 

Noise  Figure 

Receiver 

4.1  dB  (max) 

Sensitivity 

-108  dBm 

A/D  Word  Slse 

12  bit 

Elevation  Beamwidth 

Antenna 

4.8°  (min) 

Asimuth  Beamwidth 

1.4° 

Power  Gain 

34  dB 

Rotation  Rate 

12,5  RPM 

4 


77540-3 


DETECTION  REGION 

Figure  1-3.  Principal  micruburst  detection  area  for  a  hypothetical  airport. 
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In  this  report)  we  discuss  signal- to-interfereno;  considerations  for  detection  of  LAWS  with  air¬ 
port  surveillance  radars.  Receiver  sensitivity,  the  impact  of  partial  filling  of  the  fan-shaped  eleva¬ 
tion  beam,  Interference  from  precipitation  echoes  above  the  low-altitude  outflow  and  ground  clutter 
will  be  treated  to  provide  a  sense  for  the  capabilities  and  limitations  of  ASRs  in  this  role.  To  be 
specific,  we  will  use  parameters  of  the  ASR-9  in  all  calculations,  although  much  of  the  discussion  is 
also  applicable  to  the  ASR-8  vrhlch  has  similar  antenna  and  transmitter/receiver  characteristics. 
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n.  SIGNAL-TO- NOISE  FOH  WIND  SHEAR  DETECTION 


A.  Radar  Reflectivity  Factor  in  Thunderstorm  Outflows 

Figure  II-l,  reproduced  from  Wilson  et  a l.  (10),  shows  the  spread  in  the  maximum  radar 
reflectivity  factor  of  microbursts  detected  in  the  Denver,  CO  area  during  the  Joint  Airport  Weather 
Surveillance  (JAWS)  program.  The  median  reflectivity  factor  was  45  dBZ  and  in  10%  of  theee 
microburiti,  peak  reflectivity  was  20  dBZ  or  less.  Preliminary  analysis  of  FAA/Lineoln  Labora¬ 
tory  Observational  Weather  System  (FLOWS)  data  fron  Memphis,  TN  indicates  that  white  surtace 
radar  reflectivities  in  the  moist  sub-cloud  environment  at  Memphis  are  typically  larger  than  at 
Denver,  significant  wind  shear  may  still  occur  in  the  absence  of  intense  precipitation.  Of  the  27 
microbursts  identified  in  1984  with  the  MESONETf  network  of  surface  stations,  2  were  not  accom¬ 
panied  by  measurable  precipitation  (11).  Figure  II-2,  reproduced  from  Rinehart  and  Isaminger 
(12),  shows  the  distribution  of  maximum  low-altitude  reflectivity  for  FLOWS  microbursts  in  1985. 
Here,  the  median  reflectivity  factor  Is  55  dBZ  and  the  tenth  percentile  is  45  dBZ.  Two  of  these 
microbursts  had  maximum  surface  reflectivity  less  than  or  equal  to  20  dBZ. 

The  above  values  represent  the  maximum  surface  radar  reflectivity  in  thunderstorm  microbursts 
which  may  be  significantly  higher  than  the  reflectivity  in  the  region  of  greatest  wind  shear.  For 
example  microbursts  at  Denver  and  Memphis  respectively,  Tables  II- 1  and  II-2  list  this  latter 
reflectivity  factor  as  well  as  other  parameters  that  affect  microburst  detectability.  Table  II-l  was 
generated  by  S.D.  Campbell  (personal  communication)  using  radar  data  provided  by  the  National 
Center  for  Atmospheric  Research.  Table  II-2  was  compiled  by  M.W.  Wolfson  (personal  communi¬ 
cation)  using  FLOWS  radar  observations  during  the  1685  data  collection  season.  The  tabulated 
parameters  are:  (i)  the  magnitude  of  the  largest  wind  (approaching  or  receding);  (il)  the  depth  of 
the  outflow  over  the  volume  of  largest  wind;  (ill)  the  average  of  the  reflectivity  factors  in  the 
volumes  of  maximum  outflow  approaching  and  receding  velocity;  ard  (iv)  the  near-surface 
reflectivity  factor  in  the  downdraft  core  of  the  microburst. 


Table  II- Is  Example  Denver  Microburst  Features 


Date/ Time  (MDT) 

Radar 

Peak  Wind 
(««*') 

E3S9MI 

Outflow 

Depth  fkmj 

29  June  1982 

1606 

CP-2 

12 

0.8 

1.0 

23 

29  June  1982 

1606 

CP-4 

10 

0,6 

0.9 

11 

14 

8  July  1983 

1351 

CP-2 

8 

0.3 

0.5 

25 

23 

8  July  1083 

1351 

CP-4 

9 

0.3 

0.4 

10 

23 

13  July  1982 

1506 

CP-4 

9 

0.4 

0.4 

15  ^ 

28 

14  July  1982 

1444 

9 

0.8 

0.7 

- - 

13 

14 

t  A  network  of  meteorological  surface  lUtioni  that  tample  temperature,  relative  humidity,  atmoeplieiic  pree- 
sure,  rainfall  accumulation!  and  wind  epeed  and  direction  at  30  i  intervale. 
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REFLECTIVITY  (dBz) 


Figure  U-1. 


Correlogram  of  maximum  radar  reflectivity  for  JAWS  inlcrobursts  versus  maximum  velo¬ 
city  differential.  (Reproduced  from  Wilson  ef  of.  (10)). 


NUMBER  OF  MICROBURSTS 


Table  II-81  Summary  of  Mieroburat  Feature*  for  FLOWS  Memphis  Data  In  1086 

- 1  1 *— r-  " — -  "■7.Ti;.,'X,paBas3s: 


Date/  Time  (CST) 

Radar 

Peak  Wind 
(«•-') 

Height 

Resolution 

(km) 

Outflow 

Depth  (km) 

<LBzt0" 

30  April  1086 

1462 

UND 

8 

0.8 

0.8 

20-35 

36-40 

30  April  1086 

1544 

UND 

24 

0.3 

30 

45 

30  April  1985 

1608 

UND 

13 

0.4 

0.5 

25 

40 

26  June  1086 

1238 

FL-2 

r'-* 

0.1 

0.4 

50 

50 

20  Juns  1985 

1348 

FL-2 

12 

0.3 

0.4 

50 

45 

30  June  1985 

1705 

15 

0.5 

0.8 

HI 

50 

15  July  1985 

1357 

FL-2 

12 

0.1 

0.5 

25 

55 

23  July  1985 

1418 

FL-2 

12 

0.4 

20 

20-50 

10  August  1986 

1400 

FL-2 

14 

mm 

20-45 

25-45 

15  August  1986 

1519 

FL-2 

12 

0.1 

0.3 

20-45 

20-45 

24  August  1985 

1418 

FL-2 

12 

0.8 

0.8 

20 

40 

7  September  1985 
1619 

FL-2 

8 

0.2 

0.4 

10-35 

20-40 

8  September  1985 
1755 

UND 

13 

0.3 

0.7 

Particularly  In  the  Memphii  outflows,  radar  reflectivity  in  the  head  wind  and  tail  wind  eections 
of  the  mlcrobursts  may  be  significantly  lower  than  in  the  precipitation  core  that  generates  the 
downdraft.  Ae  an  example)  reflectivity  in  the  Memphis  microburet  on  30  June  varied  from  50  dBZ 
In  the  precipitation  core  to  5  dBZ  or  less  in  the  region  of  maximum  radial  velocity.  In  those 
Denver  mlcrobursts  considered  In  Table  11-1,  evaporation  in  the  dry  sub-cloud  layer  reduced  radar 
reflectivity  even  in  the  downdraft  core  to  values  less  than  30  dBZ. 

Low  radar  cross-sections  may  also  be  presented  by  gust  fronts,  once  these  have  propagated 
away  from  the  precipitation  that  generates  the  outflow.  Table  11-3,  condensed  from  Klingle  (13), 
gives  parameters  of  gust  fronts  measured  at  the  National  Severe  Storms  Laboratory  (NSSL)  in  Nor¬ 
man,  Oklahoma.  Minimum  radar  reflectivities  along  the  gust  front  range  from  approximately  0  to 
30  dBZ,  in  conformance  with  results  published  in  Doviak  and  Zralc  (14). 

For  comparison,  Figure  II-3  plots  the  minimum  detectable  signal  versus  range  for  an  ASR-0, 
plotted  in  units  of  weather  reflectivity.  The  computation,  described  in  Appendix  A,  assumed  that 
the  antenna  beam  is  filled  with  constant  reflectivity  precipitation.  The  dashed  curve  assumes  a 
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MINIMUM  DETECT! 


Figure  II-8. 


Minimum  detectable  weather  reflectivity  factor  versus  range  Tor  an  A8R-9.  The  three 
curve*  are  for  different  STC  function*  a*  described  in  the  text. 


Table  II-81  Example  NS8L  Gust  Front  Fcaturea 


Date/ Time  (CSTJ 

Radar 

Peak  Wind 
(mr'1) 

Height 

Resolution 

(km) 

Outflow 
Depth  (km) 

dB*a  |n 

30  April  1978 

2108t 

28 

0.4 

3.2 

38 

H 

0.2 

3.9 

33 

H 

02  May  1978 

1740 

NRO 

16 

0.3 

0.9 

50 

25 

19  June  1980 

2236 

NRO 

23 

0.4 

1.1 

32 

H 

15  May  1982 

1952 

0.9 

47 

-i 

40 

0.5 

2.8 

40 

9 

17  May  1983 

2228 

NRO 

43 

■QHI 

2.3 

43 

29 

26  April  1984 

2029f 

■  ■ 

3,6 

54 

9 

26  April  1984 

2148f 

36 

HH 

1.6 

37 

■ 

sensitivity  time  control  (STC)  function  that  varies  as  (R/RB){  where  R0  =  30  km.  This  is 
representative  of  STC  functions  used  for  aircraft  detection  at  sites  with  moderate  to  severe  ground 
clutter.  Clearly,  at  ranges  less  than  about  10  km  the  attenuation  of  this  STC  function  is  too  great 
for  detection  of  low  reflectivity  wind  shear  events. 

With  the  STC  disabled  (dotted  curve)  the  minimum  detectable  signal  is  below  0  dBZ  out  to  23 
km.  In  this  mode,  however,  receiver  clipping  from  ground  clutter  or  from  intense  weather  echoes 
would  severely  degrade  radar  performance.  Assuming  that  the  12-bit  A/D  converters  are  the  limit¬ 
ing  factor  for  receiver  dynamic  range,  maximum  uneaturated  signal  power  would  lie  along  a  curve 
displaced  upwards  no  more  than  85  dB. 

As  a  compromise  between  the  requirements  for  adequate  sensitivity  and  minimisation  of  clip¬ 
ping,  the  solid  line  shows  an  STC  function  that  varies  as  [It  /! f0)J  where  Rn  «  23  km.  This  func¬ 
tion  gives  a  constant  sensitivity  of  0  dBZ  out  to  23  km,  adequate  for  detection  of  most  thunder¬ 
storm  outflows.  As  will  be  shown  In  Section  IV,  ground  clutter  would  produce  saturation  in  a  rela¬ 
tively  small  percent  of  resolution  cells. 

The  specified  configuration  of  the  ASR-9  provides  a  separate  receiver  path  for  weather  signals 
only  when  circular  polarisation  is  employed.  Thus,  the  capability  to  utilise  sn  STC  function  that 
is  appropriate  for  wind-shear  detection  would  depend  on  the  air  traffic  controllers’  willingness  to 
accept  the  approximately  2,5  dB  loss  in  aircraft  target  cross-section  associated  with  the  use  of  cir¬ 
cular  polarisation.  One  alternative  would  be  to  provide  a  second  receiver’ for  LAWS  processing  at 
short  range.  This  would  receive  the  low  beam  signal  over  the  range  Interval  where  the  target  chan¬ 
nel  employs  the  hlgn  beam.  Thus  a  less  attenuating  STC  function  could  be  used  in  the  critical 
region  for  LAWS  detection  without  affecting  target  channel  parameters. 


t  Th«  itormi  oo  M  April  1871  sod  M  April  19M  produced  two  dirtied  gust  front*, 
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D.  Vertical  Structure  of  LAWS  Phenomena 


Observations  with  pencil  beam  Doppler  weather  radars  have  shown  that  the  divergent  outflow 
in  microbursts  is  often  shallow.  Of  those  JAWS  microbursls  where  there  were  sufficient  data  to 
identify  the  top  of  the  outflow,  00%  were  between  300  m  and  1200  m  in  depth  (10).  Radar  data 
from  FLOWS  in  1085  and  1988  Indicate  that  microburst  outflows  In  the  southeastern  United 
States  are  likewise  confined  to  a  shallow  layer;  as  indicated  In  Table  11-2  outflow  depths  measured 
near  Memphis  varied  from  200  m  to  800  m  in  the  region  of  largest  wind  speed.  Figure  11-4  plots 
example  vertical  profiles  of  radial  velocity  from  RH1  scans  of  microbursts  near  Memphis  and 
Huntsville,  AL. 

The  convergent  radial  velocity  signature  produced  by  gust  fronts  often  extends  considerably 
higher  In  altitude.  Kllngle  (13)  and  Doviak  and  Zrnic  (14)  Indicate  that  gust  frontB  are  typically 
observable  to  1000  m  -  3000  m  AGL. 


An  airport  surveillance  radar  will  suffer  a  significant  loss  in  effective  sensitivity  against  shallow 
thunderstorm  outflows,  owing  to  partial  Ailing  of  the  antenna  beam.  The  loss  In  sensitivity!  ,  rela¬ 
tive  to  that  plotted  in  figure  II-3,  is  given  by: 


BL  m 


r.  *"*’  * 

jJnr(»)DR(D)d<> 


U) 


where  Br(0) ,  BK(l>)  are  the  transmit  and  receive  elevation  antenna  patterns  and  1»  the  angle 
subtended  by  the  top  of  the  thunderstorm  outflow.  Figure  11-5  plots  the  one-way  elevation  pat¬ 
terns  for  the  ASR-O’s  low  (solid)  and  high  (dotted)  antenna  beams.  An  antenna  tilt  of  2.0"  Is 
assumed,  The  dashed  curve  is  the  "effective"  high-beam  pattern  given  that  signal  transmission  Is 
on  the  low  beam.  Figure  II-Q  plots  the  associated  beamillllng  loss  as  functions  of  range  for 
assumed  outflow  depths  of  300  m,  500  in,  1000  m  and  3000  m.  The  solid  curves  are  for  the  low 
receiving  beam  and  the  dashed  curves  are  for  the  high  beam.  Figure  II-7  combines  the  losses  for  a 
300  m  and  500  m  thick  outflow  layer  with  the  sensitivity  calculation  in  Figure  II-3.  (  An 
(/l  /23 '  n  )J  STG  function  Is  assumed.) 


Assuming  that  the  range  at  which  the  switch  from  the  high  to  the  low  receiving  beam  occurs  Is 
no  greater  than  30  km,  the  curves  indicate  that  echoes  from  a  3  kin  thick  gust  front  would  be 
reduced  by  no  more  than  4  dll.  For  a  shallow  mlrrohurst  however,  the  loss  eould  exceed  20  dU  In 
the  high  .team  at  operationally  significant  ranges.  The  plots  indicate  that  It  would  be  desirable  to 
use  the  low  receiving  beam,  oven  at  short  range,  fur  detection  of  shallow  thunderstorm  outflows, 
Fur  example,  at  20  km  range  the  minimum  detectable  reflectivity  factor  for  a  microburst  300  m  in 
depth  would  lie  12  tllJZ  In  the  low  beam  and  23  dUZ  in  the  high  beam.  In  ten  percent  of  the 
microbursts  measured  during  the  JAWS  program  in  a  dry  sub-cloud  environment,  peak  surface 
reflectivity  was  less  than  tho  Utter  value  (Figure  11-1).  As  indicated  In  Table  11-2,  reflectivity  fac¬ 
tors  In  the  region  of  greatest  wind  shear  are  sometimes  lower  than  23  dUZ,  even  when  the  sub- 
cloud  layer  Is  moist  as  in  Memphis. 


t  This  loss  applies  to  the  weather  echo  ■pectrinn  component  that  ii  returned  from  scatterers  In  the  outflow. 
The  Impscl  of  the  epectrum  component  returned  from  prncipltsllon  sciittereri  shove  the  outflow  Ii  treated  In 
flection  III 
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Figure  11-4.  Vertical  distribution  of  the  windspoeds  of  three  FLOWS  microbursts,  The  horltontai  dis¬ 
tance  from  the  radar  to  each  profile  is  indicated  by  the  associated  vertical  line,  This  lino 
also  marks  the  sero  point  for  the  radial  velocity  scale. 
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POWER  ((IB) 


One-way  elevation  antenna  pattern*  lor  an  A9U-U.  An  antenna  tilt  of  2.0"  la  assumed 
The  solid  line  is  the  pattern  for  the  active  (transmit  and  receive)  low  beam.  The  dotted 
line  is  the  pattern  for  the  receivc-only  high  beam.  The  dashed  line  Is  the  effective  high 
beam  pattern  (the  square  root  of  the  product  of  the  low  and  high  beam  patterns),  given 
that  signal  transmission  is  on  the  low  beam. 


BEAM  FI  LUNG  LOSS  (dB) 


5 


'igure  Il-fl.  Beamfllling  lose  versus  range  for  an  ASR-9.  The  different  curves  are  for  echoes  extending 
from  the  surface  to  the  indicated  heights 
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Figure  II-7.  ASR-9  minimum  detectable  weather  reflectivity  factor  including  the  beamlilling  losses  for  a 
300  or  600  m  deep  thunderstorm  outflow  layer.  An  (/i/23tm)*  STC  function  is  assumed. 


m.  INTERFERENCE  FROM  PRECIPITATION  ECHOES  ABOVE  THE  OUTFLOW 


A.  Single- Be  am  Measurement 

Owing  to  its  fan-shaped  elevation  beam  pattern,  an  airport  surveillance  radar  will  receive 
significant  scattered  power  from  precipitation  above  shallow  thunderstorm  outflows,  even  at  short 
ranges.  For  example,  a  500  m  deep  outflow  will  fill  the  radar’s  low-beam  3  dB  limits  only  to  5.7 
km  in  range.  At  greater  ranges,  the  velocity  spectrum  of  the  precipitation  echo  will  be  significantly 
modified  by  higher  altitude  scatterers  with  a  different  radial  velocity  (typically  lower  in  magnitude 
and/or  of  opposite  sign).  Mean  velocity  estimators,  for  example  the  pulse-pair  method,  will  there¬ 
fore  underestimate  the  true  outflow  velocity  when  the  power  contribution  from  these  scatterers  is 
significant  (e.g.  -5  dB  relative  to  power  from  scatterers  in  the  outflow  layer).  Full  spectral  estima¬ 
tion  (conventional  or  nonlinear)  might  resolve  the  outflow  velocity  component  at  lower  signal  to 
interference  ratios,  but  an  algorithm  for  identification  of  the  spectral  component  of  interest  is 
required.  In  any  case,  this  "clutter”  --  closely  related  to  the  beamfilling  loss  discussed  previously  — 
further  limits  an  airport  surveillance  radar’s  detection  range  for  LAWS  events. 

The  ratio  of  power  from  scatterers  in  the  outflow  region,  P(S),  to  power  from  scatterers  at 
higher  altitudes,  P(C),  is  given  by: 

P(S)  Jo,,Z{9)BT(0)BR(9)d9 

P(C)  f  (2) 

J,  Z(9)BT(e)BR(6)dO 

0 

=  (l/BL-1)-1 


where  the  beamfilling  loss,  BL,  was  defined  in  equation  (1).  Here  Z[6)  is  the  cloud  reflectivity 
profile  versus  elevation  angle.  The  second  equality  holds  when  precipitation  reflectivity  is  constant 
over  the  elevation  angle  interval  where  antenna  gain  is  significant  (  <  40°).  Obviously,  this  ratio 
decreases  as  the  beamfilling  loss  increases  (  i.e.  BL  becomes  smaller).  Thus  loss  of  sensitivity  for 
shallow  wind  shear  events,  owing  to  partial  beamfilling,  is  compounded  by  an  increase  in  the 
interference  background.  Figure  III-l  plots  the  power  ratio  of  equation  (2)  as  a  function  of  range 
for  the  outflow  depths  considered  previously  (300  m,  500  m,  1000  m  ?nd  3000  m).  The  calculation 
assumes  constant  radar  reflectivity  from  the  surface  to  iOOGG  m.  If  the  high  receiving  beam  is 
used,  this  power  ratio  becomes  less  than  unity  at  ranges  of  3,5,10  and  29  km  respectively  for  the 
above  outflow  depths.  If  the  ground  clutter  environment  permits  data  from  the  low  receiving 
beam  to  be  used,  these  ranges  increase  to  6,10,19  and  55  km. 

1.  Simulation  of  ASR  Velocity  Reports  Using  Simplified  Outflow  Models 

To  determine  an  ASR’s  detection  range  using  a  specific  velocity  estimator  -  the  pulse-pair 
method  -  we  simulated  time-series  weather  signals  corresponding  to  a  simple  model  of  the  vertical 
distribution  of  precipitation  reflectivity  and  radial  velocity  in  a  thunderstorm  outflow  and  the  ASR 
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Figure  1H-1.  Decibel  ratio  of  echo  power  from  precipitation  In  a  thundoratorm  outflow  layer  to  power 
scattered  from  precipitation  above  the  outflow,  The  assumed  depth  of  the  outflow  layer  la 
indicated  on  each  curve.  The  solid  curves  are  for  the  low  receiving  beam  and  the  dashed 
curves  are  for  the  high  beam,  The  computation  assumes  constant  radar  reflectivity  from 
the  surface  to  10  km  AGL, 


elevation  beam  patterns.  The  power  spectrum  of  the  precipitation  echo  sensed  by  a  fan-beamed 
radar  is  related  to  the  elevation  angle  resolved  power  spectra,  S(  V,f)  byt 


JarS(Vl9)BT(6)BH(»)d» 

r 

lorBT(t)Bx(»)dt 


(3) 


Convolution  with  the  Gauieian  ican  modulation  function,  •S'jw(V'),  account*  for  the  increaee  in 
spectrum  width  caused  by  the  rapidly  scanning  ASR  antenna.  A  white-noise  spectral  component 
corresponding  to  the  nobs  level  for  an  ASR-9  was  added  to  the  spectrum  of  equation  (3).  In  calcu¬ 
lating  the  nolie  level,  we  aesumed  use  of  an  {P  /23 km  )9  STC  function  as  discussed  in  section  d. 
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From  the  calculated  spectrum,  1  and  Q  time-series  weather  data  were  simulated  using  a  tech¬ 
nique  analogous  to  that  described  by  Slrmans  and  Bumgarner  (15).  This  Monte  Carlo  approach 
simulates  statistical  realisations  of  the  mean  weather  spectrum  where  the  amplitude  distribution  of 
each  frequency  component  is  Rayleigh  distributed  and  the  phase  is  uniformly  distributed  between  0 
and  2ir.  Inverse  Fourier  transformation  of  this  complex  spectrum  produces  a  time-series  with 
appropriate  statistical  attributes,  conforming  to  the  input  power  spectrum  shape.  We  used  128 
point  transforms  and  extracted  from  the  resulting  time-series  one  8-point  segment,  corresponding 
to  the  low-PRF  coherent  processing  Interval  (CPI)  for  an  ASR-9.  Pulse  pair  velocity  estimates) 

Tw''}  (<) 

were  formed  where  X  is  the  radar  wavelength,  ra  the  pulse  repetition  period,  the  s,  are  the  time 
samples  and  N  is  the  number  of  pulses  in  a  CPI. 

The  standard  deviation  of  the  resulting  velocity  estimates  is  large  owing  to  the  short  CPIs 
available  for  averaging  of  autocorrelation  lag  estimates  and  the  large  spectrum  width  associated 
with  the  vertical  shear  in  radial  velocity.  In  practice,  velocity  estimate  variance  could  be  substan¬ 
tially  reduced  by  incoherent  averaging  of  estimates  obtained  on  at  least  six  successive  antenna 
scans,  since  the  4.8  «  scan  period  is  short  in  comparison  to  the  time  scale  for  significant  evolution 
of  a  microburst  or  gust  front.  We  havo  assumed  Incoherent  averaging  of  six  independent  velocity 
estimates  in  the  simulations  presented. 

The  mean  and  standard  deviation  of  128  such  velocity  estimates  were  computed  at  1  km  incre¬ 
ments  from  1-30  km  to  illustrate  the  LAWS  detection  performance  versue  range  for  an  ASR.  Fig¬ 
ures  IU-2  and  III-3  plot  the  +/•  one  standard  deviation  limits  for  the  velocity  estimates,  assuming: 
(i)  outflow  depths  of  300  m,  500  m,  1000  m  and  3000  m  as  above; 

(11)  constant  reflectivity  of  30  dBZ  from  the  surface  to  10,000  m; 

(ill)  a  radial  velocity  of  -16  ms'1  in  the  outflow  layer,  changing  discontinuously  to  +5  ms*1  above 
this  layer. 

The  weather  spectrum  width  in  each  of  the  two  layers  is  taken  as  equal  to  a  relatively  high  value 
of  3  ms'1  to  simulate  the  effect  oi  additional  vertical  shear  In  radial  velocity  within  the  layers  con¬ 
sidered.  The  two  figures  treat  separately  the  low  and  the  high  receiving  beam.  Note  that  the  max¬ 
imum  range  considered  Is  well  beyond  operational  requirements  for  microbursts.  As  indicated  in 
the  introduction,  however,  detection  and  tracking  of  gust  fronts  at  relatively  long  range  from  an 
airport  may  be  useful. 

As  expected,  when  signal  power  is  dominated  by  returns  from  one  of  the  two  layers  considered, 
the  pulse-pair  velocity  estimates  accurately  reflect  radial  velocity  in  that  layer.  This  situation 
occurs  at  short  range  when  the  low-altitude  outflow  dominates  and  (for  the  shallower  outflows) 
again  at  long  range  where  echoes  from  the  overhanging  precipitation  are  the  stronger  component. 
At  intermediate  ranges  where  the  spectral  power  associated  with  the  two  components  1s  compar¬ 
able,  the  velocity  estimates  reflect  the  average  velocity  between  the  two  layers, 

We  will  arbitrarily  define  the  outflow  maximum  "detection  range"  as  the  range  at  which  the 
mean  of  the  velocity  estimates  equals  -5  ms'1,  the  midpoint  between  the  radial  velocities  In  the 
two  layers  considered.  For  the  low  beam,  this  range  is  7,  11,  and  19  km  for  outflow  depths  of  300, 
500  and  1000  m  respectively,  hi  the  high  beam  the  corresponding  ranges  are  3,  5  and  9  km.  As 
would  be  expected,  theie  detection  ranges  correspond  to  the  range  at  which  the  signal-to- 
lnterference  ratio  plotted  In  figure  III-l  drops  below  0  dB.  For  the  deeper  3000  m  outflow  model, 
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Figure  I1I-2.  Monte  Carlo  simulation  of  mean  ami  standard  deviation  of  pulse-pair  velocity  eiitimatna 
versus  range.  Each  plot  assumes:  (1)  uniform  reflectivity  from  the  surface  to  10  kmi  (11)  an 
outflow  layer  from  the  surface  to  the  Indicated  height  where  the  radial  velocity  component 
Is  -15  mi~*j  (III)  radial  velocity  of  5  mu at  all  altitudes  above  this  layer.  The  simulation 
assumed  that  the  low  receiving  beam  Is  used. 
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Figure  II1-3.  Aa  in  Figure  Iil-2,  except  that  use  of  the  high  receiving  beam  la  aaaumed, 


the  low  beam  would  detect  the  outflow  well  beyond  the  30  km  range  Limit  considered  whereas  the 
high  beam  would  detect  it  to  26  km  range. 

The  simple  outflow  models  used  here  suggest  that  the  elevation  beam  patterns  of  an  ASR  will 
limit  microburst  detection  with  straightforward  mean  velocity  estimators  to  ranges  of  about  20  km 
or  less.  Oust  fronts  should  be  detectable  at  much  longer  ranges,  given  that  they  may  extend  well 
above  1000  m  in  height  and  may  have  propagated  away  from  highly  reflective  overhanging  precipi¬ 
tation.  As  is  the  case  with  the  radar  sensitivity  issue  discussed  in  the  previous  section,  better 
detection  performance  would  be  obtained  if  the  ground  clutter  environment  allowed  for  use  of  low- 
beam  data  in  the  critical  region  within  10  km  of  an  airport. 

2.  Simulation  of  ASR  Velocity  Reports  Using  Radar  Measurements  of  Thunderstorm  Outflows 

In  this  subsection,  we  simulate  ASR  measurements  of  thunderstorm  outflows  using  radar  data 
from  FLOWS  operations  in  HunstviUe,  AL.  Two  events  are  considered  to  illustrate  microburet 
feature  variability  and  its  Impact  on  the  detection  performance  of  a  fan-beamed  ASR. 

The  method  of  simulation  is  similar  to  that  employed  in  the  previous  subsection.  RHI  measure¬ 
ments  of  the  reflectivity  factor,  radial  velocity  and  spectrum  width  through  an  outflow-producing 
thunderstorm  cell  were  used  to  calculate  the  velocity  spectrum  of  echoes  received  by  an  ASR  as  in 
equation  (3).  For  this  calculation,  we  assumed  that  the  shape  of  the  power  spectrum  contribution 
from  each  elevation  angle  In  the  input  data  was  Gaussian  with  the  measured  spectrum  width. 
From  the  calculated  spectra,  pulse-pair  velocity  estimates  were  simulated  as  before.  The  procedure 
was  repeated  for  successive  range  Intervals  In  an  RHI  tilt  to  simulate  an  ASR's  radial  velocity 
measurement  along  a  fixed  aslmuth  ‘(trough  the  outflow. 

The  outflows  treated  occurred  within  10  km  of  the  FL-2  radar  so  that  vertical  resolution  was 
175  m  or  better.  Although  the  wind  shear  in  these  events  was  not  as  intense  as  that  observed  in 
many  microbursts,  the  radial  velocity  gradient  in  both  cues  significantly  exceeded  2.5  x  10‘5  i‘l 
over  a  distance  of  several  kilometers;  this  gradient  value  was  used  In  the  CLAWS  program  as  an 
operational  threshold  for  Issuance  of  a  microburst  warning  (8). 

(a)  Case  li  22  September  1086 

Figure  III-4  shows  the  reflectivity  and  radial  velocity  fields  from  a  PPI  scan  of  a  thunderstorm 
over  the  FL-2  radar  In  Huntsville.  The  scan  is  at  0°  elevation  angla.  The  storm  generated 
moderate  rain  with  40  dBZ  echoes  extending  3000  m  above  the  surface.  Embedded  in  the  general 
outflow  from  the  rain  wu  an  ama  of  strong  divergence  centered  2-4  km  east  of  the  radar.  The 
maximum  radial  velocity  differential  here  was  13  ms'1  over  a  distance  of  2.5  km. 

The  vertical  structure  of  the  outflow  wu  measured  with  a  sequence  of  RHI  ucans  at  the  aslmuth 
angles  Indicated  by  the  overlay  on  Figure  III-4.  One  of  these  tilts  Is  displayed  In  Figure  III- 8. 
Ticks  on  the  horisontal  axis  are  at  2  km  intervals  and  ticks  on  the  vertical  axle  are  separated  by  1 
km.  The  downdraft  and  approaching  velocity  region  are  essentially  ove'  the  radar  so  that  negative 
radial  velocities  were  meuured  through  much  of  the  storm’s  depth  at  short  range.  The  receding 
portion  of  the  outflow  Is  more  typical  of  radar  measurements  of  microburst  winds.  A  shallow  layer 
of  receding  flow  waa  observed  beneath  a  much  deeper  volume  where  scattered  moved  towards  the 
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Figure  III-5. 


An  RHI  scan  of  the  reflectivity  and  radial  velocity  field  at  75"  azimuth.  Vertical  axis  ticks 
are  at  1  km  intervals  and  horizontal  axis  ticks  are  at  2  km  intervals. 
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radar.  As  illustrated  here  and  in  Figure  11-4,  the  receding  air  motion  extended  500  m  •  800  m 
above  the  surface. 

Figure  HI-8  plots  simulated  velocity  reports  from  an  ASR-9  along  the  asimuths  -where  RHI 
scans  were  performed.  The  asimuths  of  the  RHI  scans  and  the  range  extent  treated  in  the  simula¬ 
tions  are  indicated  with  overlays  on  figure  HI-4.  The  horisontal  range  ticks  at  2  km  intervals  in 
figure  III-8  are  matched  to  those  in  figure  III-5.  Shown  for  each  asimuth  are: 

(1)  "truth”  as  given  by  the  reflectivity- weighted  radial  velocity  field  averaged  over  the  altitude 
interval  0-250  m  AGL; 

(2)  the  velocity  field  estimated  using  the  pulse-pair  technique  applied  to  low-beam  data  from  an 
ASR-9; 

(3)  the  velocity  field  estimated  using  high-beam  data  from  an  ASR-9. 

We  have  simulated  a  single  realisation  of  the  velocity  estimates  assuming  incoherent  averaging 
over  six  successive  antenna  scans. 

Velocities  in  the  approaching  portion  of  the  microburst  would  have  been  accurately  measured 
with  either  receiving  beam  of  an  airport  surveillance  radar.  This  results  because  of  proximity  to 
the  radar  and  the  depth  of  the  approaching  \  elocity  region.  The  receding  flow  is  also  accurately 
measured  with  the  low  beam.  In  the  high  receiving  beam,  however,  echoes  from  the  layer  of 
approaching  precipitation  overhanging  the  receding  Bow  would  have  biased  the  velocity  measure¬ 
ments  by  as  much  as  5  ma  *.  Along  asimuths  near  90s  the  divergent  velocity  signature  would 
have  been  difficult  to  detect  using  the  high  receiving  beam. 

In  this  case,  then,  the  microburst  velocity  signature  would  have  been  accurately  measured  by  an 
airport  surveillance  radar  if  the  ground  clutter  environment  permitted  use  of  the  low  receiving 
beam.  Using  the  high  beam,  the  divergent  signature  would  have  been  measured  over  some  but  not 
all  of  the  azimuthal  extent  of  the  microburst. 


(b)  Case  2:  19  July  1986 


A  0”  elevation  angle  PPI  display  of  a  thunderstorm’s  reflectivity  and  radial  velocity  fields  is 
shown  in  Figure  HI-7 .  The  storm  occurred  slightly  to  the  east  of  the  FL-2  radar  in  Huntsville. 
Heavy  rain  produced  a  microburst  centered  8  km  from  the  radar  with  a  velocity  differential  of 
16  m«~‘  over  about  4  km. 


Eight  RHI  scans  were  performed  between  the  azimuths  indicated  by  the  overlay  on  Figure  III-7. 
One  of  these  is  displayed  in  Figure  III-8.  The  approaching  velocity  region  in  this  microburst  was 
extremely  shallow.  The  height  of  the  radial  velocity  zero-crossing  was  generally  300  m  or  less  with 
an  overhanging  volume  of  positive  radial  velocity  extending  to  1. 5-2.0  km  AGL.  The  receding  por¬ 
tion  of  the  microburst  was  considerably  deeper,  extending  about  1.0  km  AGL. 


Figure  III-9  shows  "truth”  and  simulated  ASR  velocity  reports  as  before.  The  negative  radial 
velocity  region  on  the  near  side  of  the  microburst  would  not  have  been  measured  with  either  receiv¬ 
ing  beam  owing  to  the  shallowness  of  the  outflow  layer.  In  spite  ^f  the  greater  range  from  the 
radar,  an  ASR  would  have  performed  better  in  measuring  the  deeper,  receding  microburst  winds. 
Here,  the  low  beam  velocity  estimates  would  have  been  accurate  to  within  a  few  meters  per  second. 
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Figure  III-7.  Sector  scan  at  0.0°  elevation  of  reflectivity  and  radial  velocity  fields  in  a  microburst  near 
Huntsville,  AL.  The  overlay  indicates  the  range  extent  and  azimuth  interval  between 
which  the  ASIt  velocity  reports  in  Figure  III-9  were  simulated. 
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Figure  HI-8.  RHI  scan  of  reflectivity  and  radial  velocity  at  101”  azimuth.  Vertical  axis  ticks  are  at  1 
km  intervals  and  horizontal  axis  ticks  are  at  2  km  intervals. 
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Figure  III-S.  "Truth”  and  simulated  ASR-9  velocity  reports  over  the  range-; 
in  Figure  III-7. 


An  estimate  of  the  velocity  gradient  in  this  raicroburat  using  ASR  data  would  not  have 
exceeded  1.5*  10'1  e"1.  This  is  well  below  the  previously  mentioned  CLAWS  operational  threshold. 
In  addition,  the  negative/positive  velocity  couplet  signature  of  the  microburst  would  not  have  been 
measured.  This  event,  then,  would  probably  have  been  missed  if  ASR  measurements  were  relied 
upon.  The  case  illustrates  that  extremely  •■hallow  thunderstorm  outflows  may  be  difficult  to  detect 
with  an  airport  surveillance  radar  even  when  they  occur  at  short  range. 

B.  Dual  Beam  Measurement 

One  potential  method  for  Improving  an  airport  surveillance  radar's  measurement  of  low-altitude 
wind  shear  Is  to  utilise  data  from  both  receiving  beam*  tot  (a)  "sharpen”  the  effective  beam  pat¬ 
tern;  and/or  (b)  eliminate  those  spectral  components  that  arise  front  higher  altitude  precipitation 
scatterers.  The  first  technique  implies  combination  of  signals  at  the  1  and  Q  level  u  in  a  mono¬ 
pulse  radar,  whereas  the  second  Involves  linear  combination  of  power  spectrum  (or  equivalently 
autocorrelation)  estimates  from  the  two  beams.  The  potential  efficacy  of  these  approaches  can  be 
readily  discerned  from  the  beam  patterns  in  figure  II-5.  Sensitivity  to  low  elevation  angle  precipi¬ 
tation  differs  by  14  dB  between  the  two  beams;  at  angles  greater  than  5*  however,  tha  difference  Is 
2  dB  or  less.  Thus,  "subtraction"  of  the  high  beam  signal  from  the  low  beam  signal  should  pro¬ 
vide  a  higher  degree  of  cancellation  at  large  elevation  angles  than  at  small  elevation  angles;  the 
effect  is  to  Increase  the  power  ratio  in  equation  (2). 

In  the  remainder  of  this  section  we  describe  and  evaluate  two  specific  techniques  for  combining 
signals  from  the  high  and  low  receiving  beams  to  extend  an  airport  surveillance  radar's  detection 
range  for  shallow  LAWS  events.  While  not  exhaustive,  these  implementations  illustrate  processing 
considerations  and  th«  magnitude  of  Improvement  that  could  be  obtained. 

As  outlined  In  figure  111-10,  both  approaches  presume  that  ground  clutter  haa  been  separately 
filtered  from  the  two  receiving  channels  prior  to  combination,  in  figure  III- 10(a),  I  and  Q  signals 
art  combined  so  as  to  maximise  the  power  ratio  of  equation  (2).  A  pulse-pair  mean  velocity  esti¬ 
mate  is  formed  from  the  resulting  signal.  In  figure  Iil-lO(b)  a  second  "orthogonal"  signal  channel 
Is  synthesised  that  minimises  the  power  ratio  in  equation  (2).  Autocorrelation  lag  eatimates  are 
formed  separately  for  the  two  receiving  channels  and  linearly  combined.  Since  the  autocorrelation 
function  and  power  spectrum  are  a  Fourier  transform  pair,  the  intent  is  to  cancel  that  part  of  the 
power  spectrum  that  Is  common  to  the  two  receiving  channels  (i.e.  scattering  from  high  elevation 
angle  precipitation).  A  pulse-pair  velocity  estimate  b  formed  from  the  combined  autocorrelation 
lag  estimate. 

For  both  methods,  the  range-dependent  weights  used  for  combining  signals  from  the  two  receiv¬ 
ing  beams  are  determined  from  a  model  of  the  thunderstorm  outflow,  Hers  ws  will  assume,  that  the 
outflow  extends  to  500  m  in  height,  that  precipitation  reflectivity  Is  constant  from  the  surface  to 
5000  m,  and  that  reflectivity  is  sero  at  higher  altitudes.  An  adaptive  scheme  for  combining  the 
two  channels  could  probably  Improve  performance  at  the  expense  of  increased  signal  processing 
complexity.  For  an  Initial  evaluation  however,  use  of  "static”  weights  was  felt  to  bn  adequate. 

1.  Coherent  Combination  of  I  and  Q  Signals 

For  both  of  the  above  methods,  we  require  complex  weights  for  the  combination  of  high-  and 
low-beam  1  and  Q  signals  that  maximise  the  power  ratio  of  equation  (2).  The  method  Illustrated 
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Figure  111-10.  Schematic.  of  two  method*  for  uilng  signal*  from  an  ASR’i  high  and  low  receiving  beam* 
to  Improve  wind  measurement*  In  low-altitude  thunderstorm  outflow*. 

In  figure  III- 10(b)  requires  another  sat  of  weights  that  mlnimlie*  this  ratio.  Only  the  difference  In 
the  amplitude  and  phase  response  of  the  channels  Is  significant!  without  loss  of  generality,  there¬ 
fore,  we  can  multiply  only  the  hlgh-beam  signal  by  range-dependent  complex  weights  and  sum  it 
with  the  low  beam  signal.  We  seek  therefore  extrema  In  the  ratio: 


LULL, 

P(C) 


As  before  0fl,  s  function  of  range,  Is  the  angle  subtended  by  the  top  of  the  thunderstorm  outflow.  r> 
Is  the  range-dependent  complex  weighting  coefficient  and  7,(1)  and  Fj(f)  are  the  hw  and  high 
antenna  beam  voltage  patterns.  The  magnitude  of  these  patterns  is  the  square  root  of  the  solid 
and  dotted  curves  in  figure  II-5  respectively.  We  Incorporate  the  elevation- an qle  dependent  phase 
difference  between  the  channels  into  Fa(f)  to  that  V{(i))  it  taken  at  real.  An  estimate  of  this  phase 
difference  Is  plotted  In  figure  QI-11.  This  function  was  calculated  by  Win  er  (10)  for  an  L-band 
air-route  surveillance  radar  (ARSR)  with  vertically  displaced  antenna  feeds.  To  a  reasonable 
approximation,  the  geometry  of  hit  calculation  scales  with  the  difference  in  wavelength  so  that  the 
calculated  phase  pattern  serves  at  a  rough  model  for  an  S-band  airport  surveillance  radar. 


Figure  111-11,  Phase  difference  venue  target  elevation  angle  of  signal*  In  the  high  and  low  receiving 
beanie  of  a  fan-beamed  FAA  radar.  The  phase  difference  results  from  vertical  displace¬ 
ment  of  the  antenna  feeds.  (Adapted  from  Winter  (18).) 


A a  shown  In  appendix  B,  the  extrema  of  equation  (5)  are  given  by  the  roots  of  a  quadratic  equa¬ 
tion  whose  coefficients  Involve  definite  integrals  of  thus  voltage  antenna  patterns,  weighted  by  the 
assumed  cloud  reflectivity  profile.  For  the  above  assumptions  concerning  outflow  depth  and  cloud 
reflectivity,  figure  111-12  shows  the  effective  one-way  power  beam  patterns) 

l|Vl(#){F1(#)+q7,(l)}||  (6) 

for  the  extremal  solutions,  plotted  at  ranges  of  S  km,  10  km  and  15  km.  These  are  directly  com¬ 
parable  with  the  low  and  high  receive  beam  patterns  plotted  in  figure  II-5,  The  angles  subtended 
by  the  assumed  outflow  top  (500  m)  and  the  overhanging  echo  top  (5000  m)  are  Indicated. 

Mismatch  In  the  amplitude  and  phase  patterns  of  the  input  beam  channels  severely  limits  the 
precipitation  clutter  rejection  that  can  be  obtained.  The  maximum  slgnal-to-lnterference  solution 
typically  achieves  good  suppression  over  only  part  of  the  angular  Interval  subtended  by  the 
Interfering  precipitation  echoes,  with  malnlobe  overlap  and  sldelobes  as  high  as  -8  dB  over  the 
remainder  of  this  Interval,  The  minimum  slgnal-to-intsrference  solution  places  a  notch  at  low 
elevation  angle  In  an  attempt  to  sero  power  from  the  outflow  layer.  The  upper  elevation  limit  for 
this  notch  decreases  with  Increasing  range  In  order  to  maintain  sensitivity  to  echoes  from  higher 
altitude  precipitation. 
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Figure  III-13  repents  the  provloue  calculation  (figure  III- 1 )  of  eignal-to-lnterference  ratio  venue 
range  using  the  elevation  reepome  patterns  for  there  two  weighted  combinations  of  high  and  low 
beam  signals.  Since  the  weights  were  computed  for  an  assumed  outflow  depth  of  500  m,  the  curves 
represent  the  achievable  extrema  only  for  this  case  Comparison  of  figure  111-13  and  figure  III- 1 
shows  that  for  outflow  depths  of  1000  m  or  less,  the  maximum  slgnal-to-interference  weighting  pro¬ 
vides  an  Improvement  of  approximately  4  dB  relative  to  use  of  the  low-beam  signal.  The  range 
within  which  the  power  ratio  in  equation  (5)  exceeds  unity  Increases  to  7.6,  13  and  26  hm  for 
outflow  depths  of  300  m,  500  m  and  1000  m.  This  represents  an  average  increase  of  3096  over  the 
range  achieved  hy  the  low  receiving  beam,  suggesting  a  comparable  increase  In  detection  range  for 
shallow  thunderstorm  outflows.  For  a  deeper  3000  m  thunderstorm  outflow,  the  synthesised  beam 
pattern  actually  decreases  the  slgnal-to-lnterference  ratio  sines  the  notch  falls  within  the  angular 
limits  subtended  by  the  LAWS  event.  There  is  still  good  slgnal-to-interference  power,  however,  so 
that  use  of  the  incorrect  model  in  computing  beam  combining  weights  should  not  degrade  the  abil¬ 
ity  to  see  deeper  outflows. 

Analogously,  the  signal-to-interfersncs  ratio  for  the  minimum  slgnal-to-interference  channel  is 
significantly  lower  than  in  the  hlgh-beam  signal,  particularly  beyond  5  km.  This  will  be  shown  to 
be  useful  for  measuring  the  precipitation  "clutter”  residue  In  order  to  subtract  this  residue  from 
the  maximum  slgnal-to-interference  channel. 

To  evaluate  these  dual-beam  outflow  velocity  estimators,  we  repeated  the  Monte-Oarlo  simula¬ 
tions  of  figures  IK-2  and  IQ-3,  substituting  the  effective  dual-beam  antenna  patterns  (equation  0)  in 
computing  the  elevation  angle  Integrated  weather  spectrum.  The  receiver  noise  level  was  raised  to 
account  for  the  weighted  combination  of  uncorrelated  noise  in  the  high  and  low  beam  receivers. 

Figure  ni-14  plots  mean  and  one  standard  deviation  limits  for  pulse-pair  velocity  estimates 
formed  from  the  maximum  slgnal-to-clutter  channel  as  in  figuit  III-10(a),  As  before  the  simulation 
assumes  constant  weather  reflectivity  over  altitude  and  a  20  ms"1  velocity  discontinuity  at  the  top 
of  an  "outflow  layer”  that  Is  varied  from  300-3000  m  in  depth.  The  previously  defined  "detection 
range”  it  now  0  km,  14  km  and  25  km  for  outUow  depths  of  300  m,  600  m  and  1000  m  respec¬ 
tively.  This  represents  an  average  Increase  of  20%  In  this  performance  measure,  relative  to  pulse- 
pair  velocity  estimates  formed  from  the  low  beam  signal.  As  lndlcetsd  above,  the  synthesised 
beam  pattern  actually  decreases  the  ratio  of  signal  to  precipitation  "clutter”  for  a  deeper  3000  m 
mesoscale  outflow,  relative  to  use  of  the  low  receiving  beam  alone.  This  accounts  for  the  larger 
outflow  velocity  estimate  bias  beyond  10  knt  (compare  figures  QI-2  and  III- 14).  The  defined  detec¬ 
tion  range,  however,  is  still  greater  than  the  30  km  range  limit  considered. 

2.  Combination  of  Autocorrelation  Lag  Estimates 

We  consider  next  the  linear  combination  of  autocorrelation  estimates  from  the  two  receiving 
channels  as  in  figure  111-  10(b).  In  Appendix  B  we  show  that  the  appropriate  weighting  1st 
£. 

-X  HMiW 

''  -  — t -  (7) 

Here,  J3|(tf)  and  B%(t)  are  the  two-way  power  beam  patterns  associated  with  the  maximum  and 
minimum  slgnal-to-interference  channels  respectively.  This  factor  is  simply  the  ratio  of  precipita¬ 
tion  "clutter”  residue  in  the  two  signal  channels. 
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Figure  111-13.  D«clb«l  ratio  of  echo  power  from  precipitation  in  a  thundsrstorm  outflow  layer  to  power 
scattered  from  precipitation  above  the  outflow.  The  aiiumed  depth  of  the  outflow  layer  ia 
indicated  on  each  curve.  The  solid  and  dashed  curves  are  for  synthesised  signals  that 
respectively  maximise  and  minimise  the  signal  to  precipitation  Interference  ratio  for  a  600 
m  deep  outflow. 
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Figure  111-14.  Monte  Carlo  almulatlon  of  mean  and  standard  deviation  of  pulse-pair  velocity  estimates 
versus  range.  The  assumptions  are  as  In  Figure  111-2,  except  that  use  of  the  synthesised 
maximum  signal  to  precipitation  Interference  channel  Is  assumed  as  in  Figure  III-10(a). 


As  previously,  figure  111- IS  plots  the  mean  and  standard  deviation  versus  range  of  velocity  esti¬ 
mates  generated  using  this  method.  The  simulations  were  again  performed  using  the  simple  thun¬ 
derstorm  outflow  model  described  above  with  assumed  outflow  depths  of  300  m,  500  m,  1000  m 
and  3000  m.  The  large  estimate  standard  deviations  result  from  the  " differencing"  of  statistically 
fluctuating  measured  quantities. 

Bias  In  the  velocity  estimates  obtained  with  this  method  result  frorni 

(a)  modeling  error  when  outflow  vertical  structure  does  not  match  .that  assumed  in  computing 
the  combining  weight  in  equation  (7)| 

(b)  statistical  fluctuation  of  the  autocorrelation  lag  estimates  for  the  two  channels  which  results 
in  Imperfect  "clutter"  cancellation  even  when  the  combining  weight  Is  appropriate. 


RANQB  (km) 


Figure  111-15.  Monte  Carlo  simulation  of  mean  and  standard  deviation  of  pulse-pair  velocity  estimate* 
versus  range.  The  assumptions  are  aa  In  Figure  111-2,  except  that  uae  of  the  two-channel 
method  ikotched  In  Figure  III- 10(b)  la  aeaumed. 


Aa  llluitrated  In  Appendix  B,  the  latter  effect  aeeme  to  dominate  for  thoae  caaea  where  the  actual 
outflow  depth  li  equal  to  or  greater  than  that  assumed  In  evaluating  equation  (7). 

Comparlaon  of  Figure  111-15  with  flgurea  111-3,  111-3  and  111-14  ihowa  that  for  outflowe  of  depth 
1000  m  or  leaa,  uae  of  the  "orthogonal”  signal  channel  results  In  outflow  velocity  estimate  bias  that 
Is  smaller  than  or  equal  to  that  achieved  with  either  high  or  tow  beams,  or  the  synthesised  max¬ 
imum  slgual-to-lnterference  channel.  The  previously  defined  "detection  range”  ie  greater  than  30 
km  for  the  300  m,  500  m  and  1000  m  deep  outflow  modele.  For  an  assumed  3000  m  outflow  depth 
however,  velocity  eetimate  bias  beyond  0  km  range  la  greater  than  that  achieved  with  even  the 
high  receiving  beam.  Analysis  of  the  source  of  this  error  (Appendix  B)  indicates  that  it  is  statisti¬ 
cal  In  nature,  resulting  from  the  "subtraction"  of  weighted  autocorrelation  lag  estimates  of  approx¬ 
imately  equal  magnitude.  The  error  could  be  reduced  by  further  averaging  of  the  single  channel 
autocorrelation  lag  estimates  before  combination,  at  the  expense  of  reduced  temporal  and/or  spa¬ 
tial  resolution. 


3.  Dual-Beam  Simulation*  with  Radar  Measurement*  of  Thunderstorm  Outflow* 


In  Figure  III-1B  and  III-17  we  simulate  ASR-9  velocity  report*  for  the  microbunt  caae  studies  of 
Section  III- A  uaing  both  of  the  dual  beam  technique*  considered  In  thi*  section.  A  flvo-polnt 
moving-average  filter  along  the  rang*  axis  wa*  used  to  reduce  the  large  estimate  variance  associ¬ 
ated  with  the  two-channel  approach  of  Figure  III-10(b). 

For  both  methods,  the  simulated  velocity  estimates  for  September  23  (figure  HI-10)  correspond 
well  to  the  pencil-beam  radar  measurement*  of  the  low-altitude  wind  field.  Here,  however,  the  low 
beam  of  an  ASR  would  also  have  accurately  measured  the  outflow  wind*  (figure  HI-0)  so  that  the 
dual-beam  technique*  do  not  provide  a  substantive  Improvement. 

Outflow  winds  In  the  shallow  microbunt  on  July  10,  however,  would  not  have  been  well  meas¬ 
ured  with  either  beam  of  an  ASR  used  separately.  Comparison  of  Figures  Ul-0  and  III-17  indicates 
that  modest  Improvement  -  relative  to  the  low  receiving  beam  -  would  result  from  use  of  the 
maximum  signal  to  interference  channel  alone.  With  the  synthesised  beam  pattern,  the  shallow 
approaching  microburst  outflow  would  havs  been  apparent  along  some  aslmuths  and  the  measured 
radial  velocity  gradient  would  have  been  Increased  by  as  much  as  70  percent.  The  resulting  esti¬ 
mate  ••  2.5*10'*  s'1  -  approaches  the  CLAWS  operational  threshold  but  remains  well  below  the 
maximum  gradient  measured  In  this  microburst  by  the  pencil  beam  FL-2  radar. 

Us*  of  the  two-channel  approach  of  figure  ill- 10(b)  would  have  resulted  In  unambiguous 
Identification  of  the  shallow  approaching  velocity  region  and  therefor*  a  clearly  recognisable  micro¬ 
burst  velocity  signature.  Not*  however,  that  while  the  sign  of  the  velocity  measurements  would 
now  indicate  the  presence  of  a  microburst,  wind  speed  in  the  approaching  region  of  the  outflow 
would  be  significantly  overestimated  In  many  gates.  Since  pilots  requirs  quantitative  Information 
on  microburst  wind  shear,  the  two-channel  measurements  simulated  her*  would  probably  be  Inade¬ 
quate  for  operational  usage. 

We  conclude  that) 

(a)  signal  combination  on  the  I,Q  level  to  synthesis*  a  more  favorable  beam  pattern  results  In 
only  modest  Improvement  In  thunderstorm  outflow  wind  measurement,  relative  to  data  from 
the  low  receiving  beam  alone.  Mismatch  In  the  amplitude  and  phase  responses  of  the  two 
receiving  beam*  limits  the  extent  to  which  Interference  from  overhanging  precipitation  can  be 
rejected.  Our  simulations  indicate  that  the  rang*  to  which  shallow  outflow  winds  could  be 
accurately  measured  would  be  Increased  by  roughly  30  perrent  using  this  method; 

(b)  us*  of  a  second  syntheelsed  signal  channel  to  measure  and  "caneel”  the  overhanging  precipi¬ 
tation  Interference  residue  might  provide  a  more  substantial  improvement  in  the  ability  to 
measure  outflow  winds.  The  simple  implementation  considered  here,  however,  would  not 
produce  consletently  reliable  measurements  owing  to  the  us*  of  static  combining  weights  In 
highly  variable  thundentorm  outflow  environments  and  to  statistical  uncertainties  in  the 
autocorrelation  lag  eetlmates  used  for  the  cancellation  procedure.  Future  work  will  investi¬ 
gate  whether  data-adaptiv*  methods  for  combining  information  from  th*  two  signal  channels 
would  provide  more  robust*  measurements  of  outflow  winds. 
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IV.  GROUND  CLUTTER 


The  desirability  of  using  low-beam  data  in  the  critical  region  near  an  airport  means  that  ground 
clutter  will  be  an  important  factor  in  limiting  the  LAWS  detection  capability  of  an  airport  surveil¬ 
lance  radar.  To  quantify  its  impact,  we  have  examined  ground  clutter  recordings  from  operational 
ASR-8s  at  Memphis  International  Airport  and  at  Stapleton  Airport  in  Denver.  Appendix  C 
describes  the  clutter  recording  system  and  the  data  reduction  procedure  used  to  compute  resolution 
cell  averaged  received  power  or,  equivalently,  the  clutter  cross-section  density  <r0.  The  resolution  of 
the  reduced  clutter  data  is  140  m  in  range  by  1.4*  in  asimuth. 

Figure  IV- 1  plots  the  clutter  returns  in  a  PPI  format.  The  data  have  been  scaled  to  an 
equivalent  weather  reflectivity  factor  using  equation  A-l  in  Appendix  A.  Clutter  in  the  high  and 
low  receiving  beams  are  displayed  separately  for  each  site.  Histograms  of  ground  clutter  intensity 
—  likewise  in  units  of  equivalent  weather  reflectivity  factor  —  are  plotted  in  figure  IV-2  for  the 
range  intervals  0-10  km,  10-20  km  and  20-30  km.  The  upper  abscissa  labels  give  the  correspond¬ 
ing  cross-section  density  scale  for  the  midpoint  of  each  range  interval.  Figure  IV-3  shows  the 
median  clutter  reflectivity  factor  (ensembled  over  asimuth)  as  a  function  of  range. 

Intense  ground  clutter  at  the  Memphis  site  extended  in  a  wedge  northwards  from  the  radar  and 
was  associated  with  gently  upsloping,  tree-covered  terrain  in  the  suburbs  of  Memphis.  The  max¬ 
imum  equivalent  weather  reflectivity  factor  exceeded  85  dBZ  in  the  low  beam  with  values  greater 
than  40  dBZ  measured  over  much  of  this  area.  In  the  low  beam,  less  than  1%  of  the  resolution 
cells  within  10  km  of  the  radar  were  noise  limited. 

For  comparison,  Figure  IV-4  sketches  the  runway  layout  at  Memphis  airport.  The  airport  has 
three  instrumented  runways,  two  oriented  north-south  and  one  oriented  east-west.  Secondary  run¬ 
ways  are  oriented  north-northeast  to  south-southwest  and  northwest  to  southeast.  Aircraft 
approaching  from  the  northern  sector  on  a  3*  glide-slope  would  drop  below  500  m  at  a  range  of 
9.5  km.  Thus,  throughout  the  last  stages  of  approach  where  low-altitude  wind  shear  is  a  real 
danger  such  aircraft  would  be  over  intense  ground  clutter.  Conversely,  aircraft  approaching  the 
east-west  runways  or  from  the  southern  sector  would  generally  be  over  weaker  ground  clutter  until 
the  final  few  kilometers. 

Three  areas  of  intense  ground  clutter  were  measured  with  the  Stapleton  (Denver)  ASR-8.  These 
were:  (i)  at  ranges  less  than  5  km  from  the  radar;  (ii)  an  area  to  the  north- north  west  at  15-25  km 
in  range;  and  (iii)  in  the  southwestern  quadrant  at  25-35  km  in  range.  The  second  area  is  an 
urban  region  slightly  elevated  with  respect  to  the  airport  and  the  third  corresponds  to  higher  hills 
to  the  west  of  Denver.  Equivalent  weather  reflectivity  factors  within  these  regions  were  again  fre¬ 
quently  greater  than  40  dBZ  with  maxima  exceeding  85  dBZ. 

As  sketched  in  Figure  IV-5,  Stapleton  airport  has  four  major  runways,  two  oriented  east-west 
and  two  oriented  10  degrees  clockwise  from  north-south.  Three  of  the  four  runways  are  equipped 
for  instrumented  landings.  The  principal  approach  directions  intersect  regions  of  intense  ground 
clutter  only  to  the  west  at  ranges  greater  than  30  km  and  within  5  km  of  the  radar.  Thus,  from 
the  viewpoint  of  low  altitude  wind  shear  detection,  vulnerable  aircraft  at  low  altitude  would  be 
over  severe  competing  ground  clutter  only  within  5  km  of  the  ASR-8. 
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Figure  IV- la.  PPI  display  of  ground  clutter  returns  from  the  ASR-8  at  Memphis  International  Airport. 

High  and  low  beam  returns  are  shown  separately.  Received  power  is  scaled  to  represent 
the  equivalent  weather  reflectivity  factor  using  a  fiiled-beam  assumption. 
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Figure  IV- lb.  PPI  display  oT  ground  clutter  returns  from  the  ASR-8  at  Stapleton  Air¬ 
port  in  Denver. 
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ASR-8-  MEMPHIS  LOW  SEAM  ASR-8:  MEMPHIS  HIGH  BEAM 

CLUTTER  CROSS  SECTION  DENSITY  <dB//1  m2/m2)  CLUTTER  CROSS  SECTION  DENSITY  <dB//1  m2/m2) 


Figure  IV-2a.  Histograms  of  ground  clutter  intensity  at  the  ASR-8  at  Memphis  International  Airport. 

Separate  histograms  are  shown  for  both  receiving  beams  and  for  each  of  three  range  inter¬ 
vals.  The  data  are  scaled  to  an  equivalent  weather  reflectivity  factor  (dBZ).  Upper 
abscissa  label  is  the  corresponding  clutter  cross-section  density,  a0,  at  the  midpoint  of  the 
range  interval  considered. 


59 


Figure  IV-3a,  Median  ground  clutter  inteneity  (emembled  over  ulmuth)  venue  range  for  the  ASR-8  at 
Memphis  International  Airport. 
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Figure  IV-3b.  Median  ground  clutter  lntenelty  (ensembled  over  aiirnuth)  vereui  range  for  the  ASR*B  at 
Stapleton  Airport. 
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At  ihort  rang*,  the  upper  dynamic  rang*  limit  for  thaaa  clutter  meaaurementa  is  greater  than 
would  apply  if  the  (R  /33km)1  STO  function  suggested  for  wind  shear  detection  in  Section  I!  was 
used.  Thus  we  can  compute  the  fraction  of  .eeolution  calls  where  ground  clutter  would  result  in 
receiver  clipping  if  the  suggested  STO  function  was  used.  As  a  function  of  range,  the  fraction  of 
asimuthal  cells  where  clipping  would  occur  is  plotted  in  Figure  IV-6.  At  Memphis,  both  receiving 
channels  would  be  clipped  at  almost  all  asimuths  within  about  S00  m  of  the  radar,  fiver,  in  the 
low  beam  however,  the  percentage  of  asimuth  ceils  that  would  be  dipped  drops  below  3096  beyond 
1  km  in  range.  No  clipping  would  occur  beyond  4  km.  At  Denver,  a  setup  error  for  the  recording 
system  resulted  In  loss  of  low»beam  data  within  3  km  of  the  radar.  From  3  km  to  4  km  less  than 
1096  of  low  beam  resolution  ceils  would  be  dipped  and  no  dipping  would  occur  at  greater  ranges. 
We  eondude  that  in  these  clutter  environments,  an  STO  function  that  provides  adequate  send* 
tlvlty  for  detection  of  wind-shear  associated  with  low  reflectivity  echoes  (5  *  30  dBZ)  could  be 
employed  in  both  receiving  beams  of  an  airport  surveillance  radar.  More  attenuation  than  con* 
sidered  would  obviously  be  required  for  the  flrat  few  range  gates  but  loss  of  sensitivity  at  short 
range  is  not  u  critical  since  beamfllllng  loss  is  negligible. 


Figure  rV*0a.  Probability  of  receiver  clipping  (ensmnbled  over  asimuth)  versus  range  for  the  A8R-8  at 
Memphis.  As  discussed  In  the  text,  an  (R  /23bm)a  STO  function  is  assumed. 
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Figure  lV-6b.  Probability  of  receiver  dipping  verau*  rango  for  the  ASR-8  at  Denver,  assuming  an 
(/f/23*w)a  STC  function. 


Aa  with  the  precipitation  "clutter"  eoniid«r«d  previously,  the  low-Doppler  spectral  component 
associated  with  ground  clutter  will  bi»s  meen  velocity  eetimetee  »w ay  from  the  thundentorm 
outflow  velocity.  To  eliminate  this  bin,  clutter  filtering  --  either  in  the  time*  or  frequency- 
domains  -  ie  required  prior  to  forming  the  velocity  estimate.  The  ehort  coherent  proceeding  inter* 
vale  (OPI)  of  an  airport  eurvelllance  radar,  however,  severely  limit  the  capability  to  achieve  a  nar¬ 
row  low-Doppler  stop-band  that  will  not  effect  weather  leming  capability. 

One  approach  to  the  problem  would  be  to  conatruct  a  bank  of  bandpass  Altera  using  all  the 
pulses  In  the  8-  or  10- pulse  OPI  at  in  the  ASR-9  target  channel)  the  weather  radial  velocity  would 
then  be  estimated  by  determining  the  filter  with  the  largest  output.  In  those  Alters  significantly 
removed  from  sero  Doppler,  ground  clutter  could  readily  be  attenuated  by  40  dB  or  more.  The 
drawbacks  to  this  approach  aret  (a)  the  filter  passbands  would  be  broad  (  approximately  8  ms'1) 
resulting  in  poor  velocity  resolution  for  weather  signals  and  a  broad  region  around  sero  Doppler 
where  weather  echoes  would  not  be  detected)  and  (b)  it  might  be  difficult  to  identify  that  filter 
with  peak  response  owing  to  the  random  nature  of  the  weather  signal.  Anderson  (17)  considers  the 
application  of  this  approach  to  wind  shear  detection  with  an  ASR. 

Slower  scanning  meteorological  Doppler  radars  frequently  apply  a  time-domain  digital  filter  to 
each  CPI,  followed  by  pulse-pair  velocity  estimation.  This  technique  does  not  carry  over  In  a 
straightforward  manner  to  an  airport  surveillance  radar  because  the  filter  Impulse  response  length 
must  be  almost  as  long  as  the  available  CPI  to  achieve  a  useful  transfer  function.  As  a  result, 
insufficient  filtered  data  are  available  for  velocity  estimation.  Anderson  (17)  describes  a 
modification  to  this  approach  where  the  filter  would  operate  across  two  or  three  adjacent  CPIs, 
thereby  providing  sufficient  time  samples  for  cluttar  filtering  and  velocity  estimation.  His  simula¬ 
tions  indicate  that  reliable  weather  radial  velocity  estimatee  can  be  obtained  whan  the  signal  to 
clutter  ratio  is  approximately  -40  dB  or  greater. 

An  alternative  approach  to  clutter  filtering  and  weather  velocity  estimation  with  an  airport  sur¬ 
veillance  radar  would  involve  full  spectral  estimation,  probably  using  parametric  or  "high  resolu¬ 
tion"  techniques.  The  use  of  such  estimators  ie  attractive,  given  the  small  number  of  data  samples 
available  for  resolving  the  weather  and  ground  clutter  spectral  components.  Computational 
requirements  of  these  algorithms  are  high,  however,  and  an  additional  computation  of  mean 
weather  radial  velocity  (for  example,  a  power  weighted  average  over  the  spectral  interval  outside 
the  ground  clutter  region)  would  be  required. 

Clearly  any  method  of  clutter  filtering  will  Involve  a  tradeoff  between  achieved  cluttar  suppres¬ 
sion  and  the  width  of  the  low-Doppler  "blind"  Interval  where  weather  echoes  will  be  subject  to 
significant  attenuation.  To  illustrate  reasonable  requirements  on  subclutter  visibility  at  the 
Memphis  and  Denver  airport  sites,  Tables  IV-1  through  IV-4  give  the  percent  of  radar  resolution 
cells  where  echoes  from  a  30  dBZ  microburet  would  exceed  the  interference  background  by  at  least 
8  dB,  usuming  that  ground  clutter  attenuation  of  0,  10,  30,  30  or  40  dB  is  achievable.  The  5  dB 
slgnal-to-lnterference  ratio  represents  a  minimum  value  for  reliable  weather  velocity  estimation 
with  a  pulse-pair  estimator.  The  percentages  are  calculated  within  2  km  range  intervals  out  to  10 
km,  and  then  in  10  km  Intervals  to  30  km.  The  upper  figure  is  for  the  high  receiving  beam  and 
the  lower  figure  la  for  the  low  beam. 
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Table  IV-li  Percent  Resolution  Cell*  with  31gn«l/(Clutter+Noiae)  >  6  dB 
Mtmphit  International  Airport 
Outflow  Reflectivity  ~  SO  dB* 


Outflow  Depth  =  900  m 


Range  (km) 

Clutter  Attenuation  (dB) 

Resampling 

Requirement 

0 

lJO 

20 

30 

40 

0.2 

1 

0 

1 

72 

37 

83 

65 

7 

2-4 

H 

34 

8 

87 

26 

94 

61 

1 

H 

21 

4-6 

44 

2P 

68 

46 

90 

68 

98 

93 

36 

6-8 

3 

1° 

40 

33 

60 

67 

96 

50 

8-10 

H 

42 

73 

84 

10-20 

0 

18 

0 

48 

0 

69 

T 

84 

0 

05 

i 

20-30 

0 

11 

0 

42 

0 

59 

0 

67 

0 

70 

* 

A  rerage  area  of  rauar  resolution  cells  is  greater  than  Cartesian  cell. 


Table  IV»8i  Percent  Resolution  Celle  with  Signal 


Mtmphit  International  Airport 
Outflow  Reflectivity  =**  SO  dB* 


Range  (km) 

Clutter  Attenu»Mcn  (dB) 

Reaampling 

Requirement 

30 

40 

0-2 

2 

0 

73 

37 

83 

65 

7 

2-4 

17 

2 

HE1 

98 

53 

100 

88 

21 

4-0 

41 

14 

65 

32 

100 

96 

36 

8  8 

52 

19 

77 

44 

90 

i  fi_e 

100 

99 

50 

|| 

SO 

23 

64 

i(K20 

59 

46 

98 

93 

100  " 

100 

* 

20-30 

46 

73 

80 

91 

_ Hi _ 

65 

98 

67 

99 

87 

100 

* 
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Average  area  of  radar  resolution  cell*  ii  greater  than  Cartesian  cell. 
Low  beam  clutter  data  not  available  inside  2  km. 


For  the  two  iite«,  Table*  IV-1  and  IV-3  ueume  beamfllling  lotiee  appropriate  to  a  thunder¬ 
storm  outflow  that  is  300  m  in  depth,  whereas  IV-2  and  IV-4  assume  an  outflow  depth  of  1000  m. 
While  receiver  noise  was  included  In  -alculatlng  ihe  interference  term,  we  did  not  Include  a  contri¬ 
bution  from  precipitation  echoes  above  the  outflow  layer. 

At  Memphis,  clutter  attenuation  of  20-40  dB  would  be  necessary  in  the  critical  region  within  10 
km  of  the  airport  if  more  than  half  of  the  resolution  cells  are  to  contain  usable  low-beam  weather 
data.  For  corresponding  fractional  visibility,  the  high  beam  signal  would  require  10-20  dB  less 
attenuation  at  ranges  less  than  8  km;  beyond  9  km  however,  brvmfillin;,  lose  for  the  300  m  deep 
outflow  considered  in  Table  IV-1  reduces  hlgh-beam  echo  power  to  the  point  that  system  noise 
would  prevent  detection.  A  corresponding  reduction  in  low  beam  signal  power  occurs  at  27  km 
range,  producing  the  smaller  fractional  visibility  values  In  the  20-30  km  range  interval. 

At  short  range,  attenuation  requirements  for  a  deeper  1000  m  outflow  at  Memphis  (Table  IV-3) 
are  similar  to  the  above.  However,  high  beam  data  would  now  be  usable  to  at  least  20  km  with 
20-30  dB  clutter  attenuation,  owing  to  much  reduced  beamfllling  toss. 

Near-in  ground  clutter  at  Denver  was  less  intense  than  at  Memphis.  Within  10  km  of  the 
radar,  the  percent  of  unobscured  resolution  cells  is  therefore  higher  than  at  Memphis,  given  the 
same  assumed  ground  clutter  attenuation.  For  low-beam  detection  of  a  20  dBZ  outflow  in  half  of 
the  resolution  cells  inside  10  km,  10  -  20  dB  ground  clutter  attenuation  would  suffice  even  If  the 
outflow  Is  shallow  (Table  IV-3).  A*  will  be  shown  subsequently,  however,  the  Intense  ground 
clutter  to  the  west  of  the  radar  at  ranges  greater  than  IS  km  would  result  in  large  areas  of  obscu¬ 
ration  unless  30-40  dB  attenuation  were  achieved. 

To  aid  in  Interpretation  of  the  tables,  the  last  column  shows  the  percentage  of  usable  resolution 
cells  required  for  resampling  the  radar  data  to  a  lower  resolution  (at  short  range)  Cartesian  grid. 
Because  the  "targets"  of  interest  hsre  are  extensive  relative  to  a  single  range-asimuth  cell,  detection 
of  thunderstorm  outflows  need  occur  only  in  a  relatively  small  fraction  of  available  reiolutlon  cells. 
This  approach  requires  implementation  of  a  ground-clutter  "map”  to  Identify  thoee  resolution  cells 
that  should  be  ueed  for  weather  processing.  We  chooee  200  m  x  200  m  as  the  maximum  Cartesian 
bin  else  that  would  permit  reliable  identification  of  a  1-2  km  radius  microburst.  For  this  bin  slae. 
the  last  column  In  the  tables  shows  the  percentage*  of  unobscured  resolution  cells  required  for  an 
average  of  one  usable  range/asimuth  gate  per  Cartesian  cell.  Comparison  with  the  table  suggests, 
for  example  that  this  use  of  "Interclutter"  visibility  might  permit  as  little  a*  20  dB  clutter  suppres¬ 
sion  to  be  employed  for  low-beam  data  inside  10  km  at  Memphis.  A*  Illustrated  subsequently, 
however,  the  non-uniform  asimuthal  distribution  of  the  ground  clutter  will  require  that  average 
fractional  visibility  be  significantly  higher  than  the  values  given  if  clutter  obscuration  Is  to  be  elim¬ 
inated. 

The  above  calculations  are  repeated  in  Tables  IV-5  through  IV-8  using  the  same  assumptions 
except  that  outflow  reflectivity  is  taken  to  be  45  dBZ.  in  this  case,  fractional  visibility 
significantly  exceeds  the  average  requirement  for  Cartesian  resampling  In  both  receiving  beams, 
even  without  cluttor  filtering. 
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71 


I  98  90 

._J _ SJL _ 


100 

100 

100 

100 

100 

10  D 

100 

100 

100 

100 

100 

100 

Average  area  of  radar  resolution  cells  is  greater  than  Cartesian  cell. 
Low  beam  clutter  data  not  available  inside  2  km. 


Finally,  In  flguro  IV-7  we  plot  th«  spatial  distribution  of  the  low-  beam  signal- to*lnt*rference 
ratio  for  the  two  sitei  assuming!  (a)  outflow  reflectivity  of  20  dBZ  and  beamfilllng  loss  correspond¬ 
ing  to  a  300  m  outflow  depth;  (b)  clutter  attenuation  of  20,  30  or  40  dB;  and  (e)  selection  of  that 
range-ailmutft  cell  in  each  200  m  x  200  m  Cartesian  bln  with  the  smallest  ground  clutter  residue 
for  weather  parameter  estimator.  Grey  areas  in  the  PPI  displays  correspond  to  regions  where  the 
slgna!-to-lnterference  ratio  is  between  0  and  S  dB;  black  designates  areas  where  the  ratio  is  below  0 
dB  or  else  clutter  data  is  missing.  (Compare  with  figure  IV-1  to  determine  regions  of  missing 
data).  Ths  color  rings  are  the  effect  of  asimuth-  independent  beamflUing  loss. 

The  intense  ground  clutter  north  of  the  radar  at  Memphis  was  quite  homogeneous  spatially; 
thus,  the  fractional  visibility  requirement  In  this  area  is  much  higher  than  the  aeimuthally  aver¬ 
aged  requirement  for  resampling  listed  in  the  Tables.  As  a  result,  clutter  attenuation  about  20  dB 
greater  than  suggested  by  Table  IV-1,  or  40  dB,  is  required  for  a  signal-to-interference  ratio  large 
enough  for  weather  velocity  estimation  in  this  area.  Similarly  at  the  Denver  site,  clutter  attenua¬ 
tion  of  40  dB  Is  required  if  large  areas  of  low-beam  obscuration  to  the  west  and  northwest  of  the 
radar  are  to  be  avoided.  For  both  sites,  beamfllling  lost  becomes  the  dominant  factor  limiting 
detection  of  a  shallow  20  dBZ  outflow  beyond  26  km. 

Clutter  filtering  must  be  set  up  to  handle  the  wont-cue  thunderstorm  outflow  assumptions  (i.e. 
low  reflectivity  and  shallow  depth).  Thus  the  above  arguments  Indicate  that  at  the  sites  con¬ 
sidered  40  dB  of  ground  clutter  attenuation,  in  combination  with  the  use  of  interclutter  visibility 
techniques,  would  be  required  for  velocity  estimation  with  both  beams  of  an  airport  surveillance 
radar.  As  suggested  by  figures  rV-1  and  IV-7,  this  requirement  for  filter  attenuation  could  be  con¬ 
siderably  relaxed  over  much  of  ths  area  considered  If  the  clutter  filters  were  chosen  on  a  range- 
asimuth  dependent  basis.  Additionally,  If  received  power  In  a  resolution  cell  was  significantly 
above  the  ground  clutter  residue  stored  in  a  clear  day  "map*,  it  would  be  appropriate  to  select  a 
wider-pusband  or  all-pus  filter  for  weather  parameter  estimation  in  that  cell.  These  techniques, 
which  are  used  in  the  ASR-2  six-level  weather  reflectivity  processor,  would  minimise  impact  of  the 
clutter  filters  on  weather  velocity  estimation  while  maintaining  clutter  residuu  below  the  required 
level. 
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ATTENUATION:  20  dB 
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V.  SUMMARY  AND  PLANNED  FUTURE  WORK 


A.  Summary 

We  have  examined  the  capability  for  thunderstorm  outflow  velocity  measurements  with  airport 
surveillance  radars  through  consideration  of  the  ratio  of  outflow  signal  power  to  interference  from 
receiver  noise,  ground  clutter  and  precipitation  above  the  low-altitude  wind  shear  layer.  We  indi¬ 
cated  that  the  sensitivity  of  an  ASR  is  sufiicient  to  detect  even  low  reflectivity  (20  dBZ)  wind  shear 
events  subject  to  important  caveats.  For  shallow  outflows  (depth  leu  than  about  500  m) 
beamfilling  lou  will  severely  degrade  high  beam  sensitivity,  even  at  ranges  less  than  10  km.  This 
loss  of  outflow  signal  power  is  compounded  by  an  increase  in  the  interference  background  caused 
by  precipitation  echoes  above  the  outflow  layer.  These  problems  were  shown  to  be  leu  severe  if 
data  from  the  low  receiving  beam  could  be  used  at  ranges  where  an  ASR  would  conventionally 
employ  the  high  beam. 

The  desirability  of  using  low-beam  data  in  the  critical  region  near  an  airport  lead  to  examina¬ 
tion  of  ground  clutter  measured  with  operational  ASR-8’s  at  Memphis,  TN  and  Denver,  CO.  We 
found  that  use  of  interclutter  visibility  techniques  would  permit  low-beam  detection  of  shallow, 
weakly  reflecting  thunderstorm  outflows  if  attenuation  of  40  dB  could  be  achieved  through  clutter 
"filtering”.  For  the  sites  treated  here,  this  requirement  could  be  significantly  relaxed  over  much  of 
the  radars’  coverage  areas  through  appropriate  processor  design.  We  did  not  consider  in  detail  the 
signal  processing  problem  of  achieving  the  required  clutter  suppression  with  the  short  coherent  pro¬ 
cessing  intervals  of  the  ASR-9.  An  interim  evaluation  of  this  problem  is  given  by  Anderson  (17). 

In  many  cases,  microburst  reflectivity  will  be  sufficiently  high  that  receiver  noise  and  ground 
clutter  are  of  secondary  importance  with  respect  to  an  ASR’s  detection  capability.  In  these  cases, 
the  echo  spectrum  component  associated  with  precipitation  above  the  outflow  layer  may  be  the  pri¬ 
mary  factor  limiting  the  range  to  which  the  wind-shear  event  could  be  detected  and  the  accuracy 
of  radial  velocity  shear  estimates.  We  indicated  that,  even  using  the  low  receiving  beam,  a  max¬ 
imum  detection  range  of  about  20  km  would  be  expected  for  microbursts  if  the  reflectivity  of 
overhanging  precipitation  is  roughly  equal  to  reflectivity  in  the  outflow  layer.  For  very  shallow 
outflows  (500  m  or  less),  this  range  could  be  as  small  as  5-10  km.  One  of  the  microburst  case  stu¬ 
dies  presented  illustrated  the  difficulty  of  detecting  a  shallow  outflow,  even  at  short  range. 

We  evaluated  the  extent  to  which  a  processor  architecture  that  combined  signals  from  the  high 
and  low  receiving  beams  could  produce  more  accurate  low-altitude  velocity  estimates.  In  simula¬ 
tions  with  simple  outflow  models,  coherent  combination  of  time-series  signals  from  the  two  beams 
to  synthesise  a  more  favorable  elevation  beam  pattern  resulted  in  a  30%  increase  in  the  maximum 
detection  range  for  "outflows”  less  than  1000  m  in  depth.  An  alternate  approach  that  includes 
cancellation  of  the  overhanging  echo  spectrum  component  using  an  "orthogonal”  signal  channel 
might  provide  a  more  significant  improvement  in  the  ability  to  measure  radial  velocity  in  shallow 
outflows.  Simulations  indicated,  however,  that  improvement  in  the  accuracy  of  radial  velocity  esti¬ 
mates  using  a  simple  implementation  of  this  technique  was  limited  by  the  use  of  static  combining 
weights  and  by  statistical  fluctuations  of  the  weather  signal. 

B.  Conclusions 

Overall,  our  analysis  indicates  that  an  airport  surveillance  radar  could  serve  as  a  useful  sensor 
for  measurement  of  low  altitude  wind  shear  in  the  critical  region  within  10  km  of  an  airport.  This 
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presumes  that  the  required  clutter  filtering  and  velocity  estimation  can  be  performed  using  the 
short  coherent  processing  interval  of  the  fast-scanning  radar.  Detection  of  thunderstorm  outflows 
and  accurate  estimation  of  radial  velocity  shear  at  ranges  greater  than  10  km  will  depend  on  the 
depth  of  the  outflow  and  the  distribution  of  cloud  reflectivity  in  the  vertical.  Clearly  the  detection 
probability  for  an  ASR  will  be  significantly  lower  than  for  an  appropriately  sited  pencil-beam 
Doppler  weather  radar  such  as  TDWR.  As  indicated  in  the  introduction,  an  appropriate  role  for 
airport  surveillance  radars  is  as  adjunct  wind-shear  sensors,  providing  for  example  an  alternate 
aspect  angle  to  near-  or  on-airport  outflows.  The  operational 'use  of  wind-shear  data  from  ASRs 
sited  at  airports  that  do  not  have  a  TDWR  must  be  based  on  a  reliable  estimate  of  detection  and 
false-alarm  probabilities  for  the  site  under  consideration. 

C.  Near  Term  Plans  for  Further  Analysis 

Additional  analysis  of  pencil  beam  weather  radar  data  will  provide  more  extensive  statistics  on 
the  vertical  structure  of  microbursts  and  gust  fronts.  Data  from  the  FLOWS/MIST  experiment  in 
Huntsville,  AL  and  from  the  JAWS  project  near  Denver  are  being  used  to  characterise  the  outflow 
reflectivity  distribution  in  "wet”  and  "dry”  microburst  environments  and  the  impact  of  overhang¬ 
ing  precipitation.  RHI  scans  of  microbursts  within  10-20  km  range  provide  the  necessary  resolu¬ 
tion  in  altitude. 

Volume  scans  with  sector  tilts  at  contiguous  elevation-angles,  or  RHI  scans  at  adjacent  azimuth 
angles  will  be  used  to  simulate  ASR  processing  over  an  entire  wind-shear  event,  as  opposed  to  the 
selected  azimuthal  "cuts”  considered  in  this  report.  When  dual-  or  muItiple-Doppler  coverage  iB 
available,  three-component  wind  vector  analysis  would  allow  the  data  field  to  be  translated  to 
arbitrary  ranges  and  aspect  angles  with  respect  to  an  "ASR”.  Our  intention  is  to  generate  a 
number  of  microburst  ana  gust  front  cases  for  use  in  developing  and  evaluating  signal-processing 
and  wind-shear  detection  algorithms. 

Future  simulations  will  include  the  effect  of  ground  clutter  and  the  necessary  Doppter  filtering. 
I  and  Q  data  from  the  ASR-8’s  at  Memphis  and  Denver  will  be  used  to  provide  clutter  signals  with 
realistic  amplitudes  and  temporal  statistics.  Our  first-cut  scheme  for  clutter  filtering/velocity  esti¬ 
mation  will  involve: 

(a)  interpolation  of  the  block-staggered  PRF  signal  samples  (weather  plus  clutter)  to  a  uniform 
sampling  interval.  This  would  result  in  usable  coherent  processing  intervals  two  to  three 
times  longer  than  the  eight  or  ten  available  from  the  ASR-9,  provided  that  the  weather 
echoes  are  not  aliased; 

(b)  clutter  filtering  with  a  weighted  Chebyshev  high-pas3  filter  designed  for  the  interpolated  sam¬ 
pling  frequency; 

(c)  pulse-pair  velocity  estimation  using  the  low  receiving  beam  or  dual-beam  techniques  as 
described  in  this  report. 

This  scheme  is  analogous  to  one  of  those  described  by  Anderson  (17)  except  that  it  effectively 
applies  a  shift-variant  filter  to  the  original  block-staggered  data. 

Owing  to  the  fast  scanning  ASR  antenna  and  corresponding  short  CPI,  clutter  filtering  based  on 
linear  processing  techniques  will  necessarily  produce  a  wide  low-Doppler  interval  where  weather 
echoes  will  undergo  significant  attenuation.  The  resulting  bias  in  velocity  estimates  could  be 
significant  for  LAWS  detection  when,  for  c  ample,  one  component  of  a  microburst  outflow  has  a 
low  absolute  radial  velocity  owing  to  overall  storm  translation.  Techniques  brsed  on  high- 
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resolution  ipectrtl  estimation  will  be  evaluated  to  determine  whether  the  necessary  clutter  attenua- 
tion  can  be  achieved  while  maintaining  sensitivity  to  low  radial  velocity  weather  echoes. 

As  indicated  in  this  report,  the  ground  clutter  distribution  at  Memphis  was  c  nsiderably  more 
severe  in  relation  to  aircraft  takeoff  and  departure  paths  than  at  Denver.  It  would  be  desirable  to 
obtain  ASR  clutter  data  from  additional  major  air  terminals  to  obtain  more  extensive  statistics  on 
clutter  reduction  requirements. 

FAA  usage  of  Doppler  weather  radar  data  requires  computer  recognition  of  low-level  wind  shear 
to  alert  controllers  to  the  hasard.  The  algorithms  under  development  for  TDWR  (1,2,3)  can  serve 
as  a  starting  point  for  a  corresponding  capability  for  airport  surveillance  radars.  However,  sub¬ 
stantial  modifications  will  be  required  for  the  ASR  owing  to  the  lack  of  vertically  resolved  data, 
the  generally  less  reliable  outflow  velocity  estimates  and  the  loss  of  weather  data  caused  by  clutter 
filtering.  In  addition,  the  on-airport  siting  of  ASRs  means  that  operationally  important  micro- 
bursts  may  sometimes  occur  over  the  radar,  thereby  involving  diametric  radial*.  Existing  algo¬ 
rithms  utilising  radial  velocity  shear  as  a  detection  parameter  will  require  modification  for  this 
cut 

To  provide  real  storm  data  for  assessing  LAWS  detection  with  an  airport  aurveillance  radar, 
Lincoln  Laboratory  has  deployed  a  modified  ASR-8  in  Huntsville,  AL.  A  colocated  pencil-beam 
Doppler  weather  radar  servee  as  "truth"  for  these  measurements.  The  transmitter  control  of  the 
ASR-8  has  been  modified  to  provide  the  block-staggered  PRF  waveform  that  will  be  used  by  the 
ASR-8.  I  and  Q  samples  from  both  receiving  beams  are  recorded  synchronouely  on  hlgh-denslty 
digital  tape  for  off-line  analysis. 

These  data  will  allow  for  confirmation  of  results  from  simulation  and  analysis  and  will  facilitate 
Identification  of  hardware  Issues.  While  the  number  of  mlerobursts  or  gust  fronts  occurring  at 
close  range  to  the  radar  will  be  limited,  we  expect  that  sufficient  cues  will  be  observed  during 
operatlone  In  1087  to  provide  at  leut  a  preliminary  demonstration  of  the  capabilities  and  limita¬ 
tions  of  ASRs  for  low-altitude  wind  shear  detection.  We  anticipate  that  the  combination  of  results 
from  this  testbed  radar  and  analysis  u  described  above  using  the  more  extensive  data  available 
from  meteorological  field  programs  will  provide  the  necessary  buis  for  determining  the  operational 
role  of  airport  surveillance  radars  u  LAWS  sensors. 
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APPENDIX  Ai  CALCULATION  OF  THE  MINIMUM  DETECTABLE  WEATHER 
REFLECTIVITY  FACTOR 


The  sensitivity  of  the  ASR-9  receiver  is  specified  as  >108  dBm.  This  figure  presumably  includes 
receiver  losses,  thermal  noise,  processing  guns  for  bandlimited  signals  against  noise  and  threshold¬ 
ing  requirements  for  detection.  For  comparison,  thermal  noise  at  293 *K  is  -110  dBm  for  the 
ASR-9  receiver.  We  use  the  more  conservative  figure  since  it  more  accurately  reflects  the  limit  for 
confident  detection  of  a  weather  echo. 


The  equation  for  received  power  from  precipitation  (for  example,  equation  4.30  in  Doviak  and 
Zrnic  (14))  is  then  inverted  to  calculate  the  minimum  detectable  reflectivity  factor: 

^(mm'/m1)  =  2.50  *  10”  (A.l) 


=  1.92  *  10 


PtG31\IStct**\Kw\2 

P3  *!_ 


l3TC 


If  the  STC  attenuation  is  constant,  this  limit  increases  with  range  from  the  radar  a a  R2.  The 
parameters  in  equation  (A.1)  are  given  below. 


X 

wavelength  (cm) 

11 

R 

range  (km) 

P» 

minimum  detectable  receiver  power  (mW) 

10(-lM/t°) 

Pt 

transmitted  power  (W) 

1.12  x.  10( 

G 

antenna  gain 

10(m/io) 

l 

loss  due  to  atmospheric  attenuation 

1.0 

Ir 

loss  due  to  finite  receiver  bandwidth 

l0(-2.3/10) 

IsTC 

attenuation  from  STC  function 

T 

radar  pulse  duration  (/*s) 

1.03 

6 

elevation  beamwidth  (deg) 

4.8 

* 

asimuthal  beamwidth  (deg) 

1.4 

\Kw\ 

magnitude  of  index  of  refraction  for  water 

0.93 
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APPENDIX  Bi  CALCULATION  OF  WEIGHTING  COEFFICIENTS  FOR  COMBINA¬ 
TION  OF  HIGH  AND  LOW  BEAM  SIGNALS 

1.  Combination  of  I  and  Q  time-series  signals 

We  seek  first  a  range-dependent  complex  weight: 

ij(fl)  =  a(R)  +  i/3{R)  (B.l) 

such  that  the  ratio  of  signal  power  from  a  low  altitude  outflow  layer  to  interference  from  echoes 
above  the  outflow  layer: 

P(S)  /a,,^(*)liv1(#){^)+nva(fl)}||ad» 

P(C )  i-  (B,2) 

4  zmVimvM+nVtmi3** 

is  an  extremum.  The  variables  in  (B.2)  were  defined  in  the  discussion  of  equation  (5)  in  the  text. 
The  elevation-angle  dependent  phase  difference  between  the  high-  and  low-beam  voltage  antenna 
patterns  is  incorporated  into  the  high-beam  pattern: 

v,(*)  =  s{9)  +  <«(#)  (B.3) 

so  that  V\(0)  is  taken  as  real. 

We  shall  employ  the  following  definitions  in  deriving  the  extrema  of  equation  (B.2): 

Ax  =  j\(9)Vf(9)d9 
Bx  =  £  Z{9)P{9)d9 
Cx  =  f‘'z(9W(0)d9 
=  j''z(9)Vl(9)6(9)d  9 
Ex  =  Cz(9)Vx{9)K(9)d9 

9 

Aa  =  J^Z(9)V?(9)d  9 

9 

B3=  Jtfz(8)?(9)d9  (B.4) 

9 

C2  =  JtTZ[6)K2{$)d9 

0 

9 

D3  =  JtrZ(9)V  l(9)S(9)d  9 
0 

9 

=  JrZ{8)Vl(8)K(9)d9 
®0  =  A2D  \  -  A\D2 
ai  —  A2(B{+Ci)  -  Ax(B2+C2) 
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<i3  —  D3E)  -  DxE3 

Z?a(/7 1+ Ci)  -  Di(B2+Ci) 

3t  —  Ei(fla+Cj)  +  E)(B  1+ C\) 

4o  —  A i£j  -  A tE i 

These  quantities  are  entirely  defined  by  the  voltage  antenna  patterns,  an  assumed  upper  boundary 
for  the  outflow  layer  and  an  assumed  profile  of  precipitation  reflectivity. 

Differentiating  equation  (B.2)  with  respect  to  a  and  /3  and  setting  the  result  equal  to  sero  gives 
a  pair  of  simultaneous  equations: 

/ a(a^)»i(«)  -  /  =  0  (B.5) 

/a(a.£)M«)  -  /  i(«i0)V«»)  =  0 

where 

/ ,M)  -  Ax  +  2a/?,  -  MEt  +  («5  +  0J)(£,  +  C,) 

/2(o,jJ)  =  A,  +  2a£>3  -  20£,  +  (aJ  +  /Ja)(B,  +  Ct) 
p,(o)  =  2D,  +  2o(fl,  +  Ci)  (B.0) 

y  a  (or )  =  2D,  ■+*  2ot[D,  4-  C^) 

+  C,) 

'V*  +  (7j) 

It  follows  that; 

(B.7) 

*.(•)  Hf)  11 

This  can  be  expanded  to  yield  a  linear  relation  between  the  real  and  imaginary  parts  of  the  weight¬ 
ing  coefficient; 

a4a  -  Oj/J+oa  =  0  (B.8) 

Substituting  this  relation  back  into  equatirm  (B.S)  and  expanding  terms  yields  a  quadratic  equa¬ 
tion  for  the  real  pert  of  the  weighting  coefficient; 


{2(as=  f  at»)^}esa  + 


{  2  fl  |  cl  4  ~  2  d  g  (1  3 


•+  2(af  +  a 4  ) — }a  + 


{a0u4  +  60(J3  i-  2<jj - *■  a(aj}  —  0 

<tj 

The  route  uf  equation  (B.9),  in  combination  with  (13.8)  define  the  desired  complex  weights. 


2.  Combination  of  autocorrelation  lag  estimate* 

We  consider  next  the  combination  of  autocorrelation  lag  estimate*  from  two  receiving  channels 
(i.e.  the  maximum  and  minimum  signal  to  interference  channels  considered  above).  We  denote  the 
'.oruplex  time-samples  for  channel  j  by  t/j  which  can  be  resolved  into  a  component  due  to  scatter¬ 
ing  from  the  thunderstorm  outflow  layer,  tu,  and  from  precipitation  above  the  outflow,  e,^: 

xi.J  =“  1.2  (B.10) 


The  autocorrelation  functions  at  a  one  sample  delay  are  them 

*»(W)“  <  *<+i >  (B.ll) 

=“  <  'i+ijt'u  >  +  <  *l+Ue*U  > 

where  cross-terms  vanish  owing  to  the  expectation  operation.  The  desired  information  on  outflow 
radial  velocity  is  contained  in  the  phase  of  the  first  terms  on  the  right.  As  a  simplifying 
assumption,  we  take  the  weather  echo  spectrum  to  be  constant  In  altitude  within  the  two  layers 
considered.  Then  the  different  elevation  angle  weightings  for  the  two  receiving  channels  do  not 
affect  the  phase  of  the  signal  and  "clutter”  components  of  the  autocorrelation  function.  This 
allows  for  cancellation  of  the  clutter  components  in  (B.ll)  through  linear  combination  of  the  auto¬ 
correlation  samples  from  the  two  receiving  channels! 

*«(»-,  ,1)  +  7*.. (r, ,2)  —  <  **+Ms*M  >  +  1<  >  (B.12) 

The  combining  weight: 


<  «i+u4*t,a  > 


(B.13) 


is  a  real  number  by  virtue  of  the  above  assumption  on  weather  variation  with  height. 


If  we  further  assume  that  the  elevation  angle  resolved  weather  echo  spectra  are  Gaussian  in 
shape,  then  the  precipitation  clutter  component  of  the  autocorrelation  function  at  one  sample  delay 
is  related  to  the  total  power  of  this  spectral  component  byi 


*..(W)  —  <  «<+i, /4V>  > 

This  allow*  the  combining  weight  to  be  written 


*..(0.1) 

“*..(0,2) 


-tfmuvmvM  +  niwm* 

w 

jfimvmvA*)  -  <f.v”,(#)>nBflt# 


(B- 14) 


(B.16) 


where  171  and  rj3  are  the  coherent  combining  weights  associated  with  the  maximum  and  minimum 
signal-to-preeipitatlon  interference  channels  respectively.  The  combining  weight  7  for  the  auto¬ 
correlation  functions  can  again  be  computed  from  knowledge  of  the  antenna  elevation  patterns  and 
an  assumed  model  for  reflectivity  versus  altitude. 
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3.  Illustration  of  error  sources  for  dual-channel  velocity  estimates 


Aa  indicated  in  section  III,  errors  in  estimating  outflow  radial  velocity  with  the  dual-channel 
method  illustrated  in  figure  III-10(b)  result  fromt 

(a)  discrepancies  between  actual  outflow  structure  and  the  model  used  in  evaluating  (B.15)| 

(b)  statistical  error  in  estimation  of  the  autocorrelation  functions  (B.ll)  which  results  in  Imper¬ 
fect  "clutter"  cancellation  even  when  the  assumed  outflow  model  is  appropriate. 

Figure  B-l  illustrates  the  error  Introduced  by  an  Inexact  model  for  thunderstorm  outflow  struc¬ 
ture.  As  in  the  text,  combining  weights  were  computed  assuming  an  outflow  depth  of  500  m  and 
uniform  reflectivity  from  the  surface  to  5000  m.  Assuming  that  the  autocorrelation  function  can 
be  exactly  determined,  the  figures  plot  estimated  velocity  venue  tangs  when) 

(a)  actual  outflow  heights  are  300  m,  500  m,  1000  m  and  3000  m| 

(b)  reflectivity  Is  constant  from  the  surface  to  10,000  mi 

(c)  outflow  radial  velocity  le  -15  ms'1  and  overhanging  precipitation  radial  velocity  is  5  ms'1. 

The  figure  is  directly  comparable  with  figure  III- 15  in  the  text. 

It  is  seen  that  the  magnitude  of  the  error  introduced  is  email  when  the  actual  outflow  depth  is 
equal  to  or  greater  than  the  assumed  depth,  500  m.  For  a  shallower  300  m  depth  outflow,  the 
inappropriate  choice  for  the  weighting  coefficient  results  in  significant  bias  in  the  velocity  estimates 
beyond  about  6  km. 

Figure  B-2  simulates  an  autocorrelation  lag  "estimate  error”  for  each  channel.  The  calculations 
are  as  in  the  preceding  figure  except  that  an  "error  vector”  has  been  added  to  each  autocorrelation 
function.  The  magnitude  of  this  error  is  uniformly  distributed  from  0  to  the  true  magnitude  of  the 
autocorrelation  function)  the  phase  Is  taken  as  uniformly  distributed  between  0  and  2ir. 

The  mean  of  the  resulting  velocity  estimates  may  differ  significantly  from  the  values  In  figure 
B-l  when  the  weighted  magnitude  of  the  autocorrelation  function  in  the  second  receiving  channel  is 
comparable  to  that  in  the  first  channel.  In  this  circumstance,  the  contribution  from  the  "error" 
components  may  substantially  alter  the  phase  of  the  vector  formed  by  combining  the  two  func¬ 
tions.  The  resulting  velocity  estimates  will  be  subject  to  large  errors  and  may  be  distributed  over 
a  significant  portion  of  the  Nyqulst  Interval.  The  effect  is  to  bias  the  mean  of  the  velocity  esti¬ 
mates  towards  the  center  of  the  Nyqulst  interval.  Observe  that  significant  differences  between  the 
mean  of  velocity  estimates  in  figures  B-l  and  B-2  coincide  with  large  estimate  standard  deviations. 
With  the  exception  of  the  300  m  "true”  outflow  depth,  this  statistical  bias  is  generally  of  greater 
magnitude  than  the  modeling  error  treated  in  figures  B-l. 

The  simple  geometrical  simulations  used  in  this  appendix  reproduce  the  qualitative  features  of 
the  time-series  signal  simulations  plotted  in  figure  111-15;  remaining  differences  are  due  to  the  sim¬ 
plifying  assumptions  used  here  for  the  distribution  of  the  autocorrelation  lag  estimate  error. 
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Dual-channel  velocity  eitlmate  veraua  range  Including  an  error  contribution  from  etatletical 
uncertainty  In  eitlmatlon  of  autocorrelation  functions. 


APPENDIX  Ci  GROUND  CLUTTER  MEASUREMENTS  AND  DATA  REDUCTION 


Bill  Drury  at  Lincoln  Lab  has  designed  an  Instrumentation  package  to  record  ground  clutter 
with  an  operational  ASR-V  or  ASR-8.  The  system  is  sketched  in  Figure  C-l.  On  each  antenna 
scan,  1  and  Q  samples  from  a  sector  that  Is  98  range  gates  by  266  pulse  transmissions  are  buffered, 
then  written  to  magnetic  tape  during  the  remainder  of  the  Kan.  The  radar  transmits  at  a  uniform 
pulse  repetition  frequency  of  1024  i"1  while  scanning  aslmuthaliy  at  a  nominal  rate  of  12.5  RPM. 
Range  gate  spacing  for  an  ASR-8  is  140  m.  Thus  the  wedge-shaped  area  recorded  -n  one  antenna 
Kan  is  18.5*  by  13.4  km.  The  starting  asimuth  and  range  for  data  collection  are  selected  with 
thumbwheel  twitches. 

In  operation,  data  from  a  given  wedge  were  recorded  on  five  successive  antenna  scans  to  allow 
for  power  averaging  or  examination  of  ecan-to-scan  fluctuations.  The  measurements  wait)  repeated 
with  new  start  asimuths  or  ranges  in  order  to  provide  full  asimuthal  coverage  from  0-30  km.  As 
is  evident  In  Figure  IV-1,  small  range-aiimuth  sectors  were  mined  owing  to  operator  error  or 
recording  aystem  malfunction. 

The  parameters  of  an  ASR-8  are  quite  similar  to  the  ASR-9  parameters  listed  In  Table  1-1.  For 
the  ground  clutter  measurements,  an  STC  function  that  varied  as  (R/18.S  Am)3  was  used.  This 
selection  results  in  constant  sensitivity  with  range  inside  18.5  km  when  the  data  are  used  to  com¬ 
pute  the  ground  clutter  croes-section  density  Figure  C-2  plots  the  dynamic  range  limits  of  the 
measurements  in  these  units.  Inside  18.6  km,  clutter  cross-section  densities  between  -72  dB 
(//  lm8/ms)  and  -2  dB  can  be  measured  If  data  from  the  two  receiving  beams  are  merged. 

Table  0-1  lists  the  sites  where  measurements  have  been  made  as  well  as  the  antenna  heights 
and  prevailing  wind  speed.  The  radar  at  the  FAA  Technical  Center  wu  used  to  test  the  equip¬ 
ment  and  to  develop  operating  procedures.  Measurements  at  Huntsville  were  made  primarily  to 
support  an  experimental  program  in  whleh  Lincoln  Laboratory  is  involved.  In  this  report,  there¬ 
fore,  we  considered  only  data  from  radars  at  the  major  air-terminals  in  Memphis  and  Denver. 


Table  C-li  Sites 

'or  Ground  Clutter  Measurements  with  ASRs 

Sit « 

Radar 

Prevailing  Wind  Speed 

FAA  Technical 
Center,  Atlantic 

City,  N.J. 

ASR-8 

<  5  m#"1 

77  feet 

Memphis  Interna¬ 
tional  Airport 

ASR-8 

<  5  ms'1 

87  feet 

Madison  County 
Airport,  Hunstvllle, 
AL 

ASR-7 

calm 

. 

57  feet 

Stapleton  Airport, 
Denver,  CO 

ASR-8 

calm 

17  feet 

Data  reduction  consisted  oft  (i)  power  averaging  of  returns  in  each  resolution  cell  over  the  five 
succeaelve  antenna  scans  and  the  20  pulses  transml.'ed  as  the  antenna  swept  across  one 
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Figure  C-l.  ASR  ground  clutter  recording  »y»tem, 
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Figure  C-2,  Dynamic  range  limits  for  ASR  ground  clutter  measurements  versus  range  from  radar. 
Expressed  In  terms  of  clutter  cross-section  density,  i*0' 


beamwldthj  (11)  sealing  of  these  averaged  sampled  intensity  values  to  represent  received  power  In 
milliwatts.  The  scaling  constant  was  calculated  by  comparing  ths  noise  level  for  the  sampled  data 
to  the  nominal  receiver  noise  level  for  an  ASR-8.  Data  from  the  furthest  4.2  km  range  Interval 
were  averaged  to  estimate  the  sampled  noise  level)  visual  examination  of  a  PP1  display  was  used  to 
exclude  sectors  within  this  interval  containing  ground  clutter  The  averaged  power  levels  in  each 
resolution  cell  were  then  thresholded  4  dB  above  this  noise  level  to  exclude  receiver  noise  from  the 
analyst. 

Thle  report  considered  therefore  the  sean-avernged  ground  clutter  Intensity.  When  a  resolution 
cell  contains  many  scatterers  of  comparable  cross-section,  relative  movement  on  the  order  of  the 
radar  wavelength  can  remit  in  large  fluctuations  in  signal  Intensity.  The  time-icale  for  such  move¬ 
ment  Is  generally  large  compared  to  the  duration  of  a  CPI  but  small  relative  to  the  scan  period. 
Clutter  Intensity  measurements  for  s  given  resolution  cell  therefore  fluctuate  on  s  scsn-to-scan 
basis  according  to  a  non-central  Chi-squared  distribution.  In  ths  limiting  case  that  no  large, 
discrete  scatterere  are  present  In  a  resolution  cell,  ths  distribution  function  is  exponential  and  the 
standard  deviation  of  the  clutter  intenelty  samples  equals  their  mean. 
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